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I  3  ABSTRACT 
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Pushing  and  incorporation  of  foreign  particles  during  solidification  were 
studied.  Trapping  was  favored  by  rapid  freezing,  a  small  amount  of  stirring, 
and  the  presence  of  gas  bubbles  at  the  interface,  and  by  formation  of  an  acute 
angle  between  the  interface  and  the  container.  With  faceted  growth  trapping 
differed  from  grain  to  grain  and  was  sometimes  favored  at  a  grain  boundary. 

The  temperature  gradient  had  no  effect  on  particle  trapping.  Particle 
mixtures  were  separated  by  slowly  increasing  the  freezing  rate.  This  is  called 
particle  chromatography.  A  theory  was  developed  for  nucleation  of  gas  bubbles 
during  crystallization. 
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ABSTRACT 

1  i 

i 

Foreign  particles  are  present  in  nearly  all  solid  organic 
chemicals,  including  laboratory-grade  chemicals.  Metals  may  similarly 
contain  particles  of  slag  and  furnace  refractories  following  extrac- 
tive  metallurgical  processes.  These  particles  are  not  only  sources 

i 

of  impurities  but  may  induce  defects  and  change  many  properties  of  a 
solid. 

I 

It  has  been  known  for  many  years  that  growing  crystals  some¬ 
times  reject  and  push  foreign  particles  in  addition  to  segregation  of 
soluble  impurities.  The  actual  interaction  between  particles  and  a 
solid-liquid  interface  was  recently  studied  experimentally  and 
theoretically,  and  some  particle  pushing  and  trapping  phenomena  ; 
clarified.  At  sufficient  low  growth  rates,  nearly  all  particles  were 
pushed  by  a  freezing  interface.  As  the  growth  rate  was  increased,  a 
critical  velocity  Vc  was  reached  beyond  which  particles  started  being 
trapped  into  the  growing  crystal.  Jhe  measured  values  of  Vc  were 
different  for  each  particle-matrix  material  system  and  possessed  no 
discernible  pattern.  Prior  to  the  present  work,  insufficient  infor¬ 
mation  was  known  about  particle  pushing  to  permit  design  of  a  separa- 

t  5 

tion  process  for  producing  particle-free  ultrapure  materials. 

Two  different  experimental  methods  were  employed  in  this  study. 
An  improved  Bridgrnan-Stockbarger  apparatus  was  first  constructed  per¬ 
mitting  microscopic  observation  of  particle  pushing  and  trapping  at 
the  freezing  interface.  Interface  shape,  growth  container,  growth 
orientation,  and  presence  of  bubbles  all  significantly  influenced 

xi  !  , 


i 


i 


For  the  carbon-naphthalene  system  Vc  was  21  mm/hr.,  but  it  was  re¬ 
duced  to  18  mm/hr.  when  the  particles  were  trapped  next  to  the  tube 
wall  with  a  convex  interface.  For  the  copper-naphthalene  system  Vc 
was  about  18  mm/hr.  The  interfacial  temperature  gradient  had  no 

discernible  effect  on  V  .  The  formation  of  bubbles  was  prevented  by 
!  C 

!*  evacuating  the  growth  container.  In  an  attempt  to  increase  Vc  by 
stirring  of  the  melt,  insertion  of  a  propeller  or  paddle  agitator  was 
'  found  to  be  impractical. 

Horizontal  zone-refining  with  rotation  was  then  adopted  to  in¬ 
crease  V  .  This  apparently  was  the  first  time  this  technique  was 
applied  for  removing  particles.  This  method  eliminated  the  detri¬ 
mental  effects  resulting  from  settling  of  the  particles  and  formation 
ov  bubbles  which  were  encountered  in  the  first  method.  Thus,  the 
measured  values  of  Vc  were  increased  300%  for  copper  and  50%  for  carbon 
particles  in  naphthalene.  This  new  separat  on  process  has  proven  to 
be  simple,  effective  and  economical  (product  cost  is  about  22  cents 
per  gram  for  particle-free  naphthalene).  Furthermore,  horizontal  zone- 
1  refining  with  rotation  and  vertical  Bridgman  technique  were  used  to 
separate  carbon  particles  from  a  mixture  of  carbon  and  copper  particles 
in  naphthalene  by  taking  advantage  of  their  difference  in  Vc>  This  new 
separation  technique  is  christened  "particle  chromatography." 

The  measured  figures  of  Vc  were  in  good  agreement  with  those 
:  calculated  by  Bolling  and  Cisse's  equation.  Theoretical  treatments  on 
formation  of  air  bubbles,  effect  of  crystal  dimension  on  the  inter¬ 
facial  temperature  gradient,  effect  of  temperature  gradient  and 
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horizontal  rotation  on  particle  pushing  were  made  in  this  study.  In 
general,  they  agreed  with  the  present  experimental  results  and 
previous  data. 

At  least  three  objectives  have  been  accomplished  in  this  dis¬ 
sertation.  (1)  Discovery  of  some  important  phenomena  leading  to 
better  understanding  of  pushing  of  particles  by  a  freezing  interface 
(2)  The  development  of  a  new  separation  process  enabling  economical 
removal  of  foreign  particles  from  organic  compounds.  (3)  The  inven¬ 
tion  of  a  new  technique  to  separate  mixtures  of  different  particles, 
opening  up  a  potentially  new  field--particle  chromatography. 
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CHAPTER  I 


INTRODUCTION 

The  demand  for  highly  pure  solid  materials  is  continuously 
growing  for  research  and  applications,  particularly  in  the  fields  of 
semiconductors  and  pharmaceuticals.  While  separation  of  soluble 
impurities  has  been  extensively  investigated  during  the  past  two 
decades,  removal  of  i nsoluble  impurities  has  been  largely  neglected. 
Foreign  particles  are  present  in  nearly  all  organic  compounds,  in¬ 
cluding  laboratory-grade  chemicals.  Metals  may  similarly  contain 
particles  of  slag  and  furnace  refractories  following  extractive  metal¬ 
lurgical  processes.  These  particles  are  not  only  source  of  impurities 
but  may  induce  defects  such  as  dislocations  [1],  and  change  many 
properties  of  a  solid. 

It  has  been  known  for  many  years  that  growing  crystals  some¬ 
times  reject  and  push  foreign  particles  in  addition  to  segregation 
of  soluble  impurities  [2].  This  observation  suggested  the  possibility 
of  removing  particles  by  solidification  processes  [3].  The  actual 
interaction  between  particles  and  a  freezing  interface  has  recently 
been  studied  [4-10],  and  some  particle  pushing  phenomena  have  been 
clarified.  At  sufficiently  low  crystal  growth  rates,  nearly  all  par¬ 
ticles  were  observed  to  be  pushed  by  the  solid-liquid  interface.  As 
the  growth  rate  was  increased,  a  critical  velocity  V£  was  reached 
beyond  which  all  particles  of  a  particular  type  were  trapped.  The 
measured  values  of  Vc  were  different  for  each  particle-matrix  material 
system,  and  possessed  no  discernible  pattern.  Prior  to  the  present 
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study  insufficient  information  was  known  about  particle  rejection  to 
permit  engineering  design  of  such  a  separation  process. 

The  purpose  ol  this  research  was  to  learn  ^  'e  about  the  funda¬ 
mental  nature  of  particle  pushing  and  trapping  phenomena,  and  to  study 
the  effects  of  some  of  the  design  variables  of  a  separation  process. 
Selected  particles  were  mixed  in  transparent  organic  melts.  In  order 
to  measure  the  critical  rate  for  trapping,  the  solidification  rate 
was  gradually  increased.  Two  different  techniques  were  employed.  An 
improved  vertical  Bridgman-Stockbarger  growth  apparatus  was  first  set 
up  for  direct  observation  of  the  particle  pushing  and  trapping  at  the 
solid-liquid  interface.  A  horizontal  zone-refiner  with  rotation  was 
then  employed.  Since  some  of  the  organic  compounds  used  here  froze 
non-faceted  like  most  metals  [11],  the  results  of  this  work  should  be 
equally  applicable  to  removal  of  particles  from  metals. 

In  the  next  chapter  the  prior  literature  is  critically  re¬ 
viewed,  and  some  theoretical  considerations  needed  for  this  work  are 
described.  This  is  followed  by  Chapter  III  which  describes  the 
detailed  experimental  methods  and  materials  employed  for  the  research. 
Chapter  IV  shows  the  experimental  results.  Chapter  V  explores 
potential  applications  with  some  experimental  evidence  and  also  gives 
economic  aspects  of  the  developed  process.  Interpretation  of 
experimental  results  is  given  in  Chapter  VI.  Accomplishments  and 
important  results  drawn  from  this  research  are  summarized  in  Chapter 
VII.  Finally,  the  extension  of  the  present  work  to  several  fields  is 
suggested  in  the  last  chapter. 


CHAPTER  II 


LITERATURE  STUDIES 

This  chapter  consists  of  two  sections:  a  review  of  previous 
work  and  some  theoretical  considerations  required  for  the  present 
work.  The  prior  literature  is  reviewed  in  chronological  order. 
Earlier  work  is  discussed  under  "Crystallization  Pressure",  since 
before  1960  many  investigators  interpreted  that  a  growing  crystal 
could  exert  a  force  which  pushed  an  external  body.  This  is  followed 
by  a  review  of  more  recent  papers  in  which  the  interaction  between 
particles  and  a  solid-liquid  interface  was  examined.  Of  these  recent 
studies,  the  work  of  Bollinr  and  Cisse  [8,9,10]  appears  most 
significant  for  particle  pushing  and  trapping.  Their  theory  is  thus 
reviewed  in  detail.  Some  theoretical  considerations  needed  for  this 
research  are  discussed  in  the  second  section  of  this  chapter,  such  as 
the  structure  of  a  solid-liquid  interface  and  the  intertacial  tempera 
ture  gradient. 

A .  Review  of  Previous  Work 
1 .  Crystallization  Pressure 

It  has  been  known  for  many  years  that  in  addition  to  segrega¬ 
tion  of  soluble  impurities,  a  growing  crystal  either  repels  or  entrap 
foreign  solids  [2,4],  Many  interesting  examples  have  been  observed 
[12].  Crystals  of  epsomite  (MgS0^'7H20)  grown  on  the  bottom  of  salt 
lakes  repelled  particles  of  mud  and  were  therefore  colorless  and 
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transparent.  On  the  other  htnc,  crystals  of  astrakanite 
(Na2S04-MgS04-4H20)  formed  under  similar  conditions  were  usually  dark 
colored  due  to  entrapment  of  mud  particles.  Gypsum  crystals 
(CaS04-2H20)  growing  in  clay  rejected  the  particles  in  their  proximity 
and  were  often  clear,  while  the  same  crystals  formed  in  sand  exhibited 
abundant  inclusions  of  sand  particles.  These  phenomena  of  particle 
repulsion  or  entrapment  may  be  related  to  a  force  exerted  by  the 
crystal  on  the  particle  during  crystalli  ation. 

Becker  and  Day  [14]  were  the  first  to  perform  experiments  in 
1905.  They  put  forward  the  view  that  a  crystal  could  exert  a  pressure 
when  growing.  They  placed  a  loaded  glass  plate  on  the  top  face  of  an 
alum  crystal.  The  crystal  continued  growing  even  when  the  load  was 
increased  to  one  kilogram.  They  attributed  the  effect  to  a  linear 
crystallization  force".  The  crystallization  pressure  was  defined  as 
the  maximum  pressure  exerted  by  a  crystal  at  which  the  growth  ceased 
at  a  given  supersaturation  {or  supercooling).  Since  then  many 
investigators  have  pursued  the  measurement  of  crystallization  pressure, 
but  have  not  agreed  nn  the  values  of  the  pressures  that  might  occur. 
Some  found  values  of  about  20  Kg/cm2  [4,14,15,16]  while  others  found 
only  a  few  g/cm2  [17,18].  The  former  used  Becker  and  Day's  method, 
the  latter  used  Shubnikov's  method  [16].  Shubnikov  used  a  growing 
crystal  to  encounter  a  glass  ball  which  was  attached  to  a  sp.ing.  The 
crystallization  pressure  was  determined  when  the  ball  began  entering 
the  crystal.  Those  who  took  the  pressure  to  be  large,  related  the 
pushing  force  to  the  energy  of  the  phase  change.  Those  who  took  the 
pressure  to  be  small  explained  the  effect  as  a  surface  interaction. 
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I  )  fact,  the  supply  of  supersaturated  solution  to  the  contact  region 
is  the  limiting  process  in  the  latter  case. 

Corre.ns  and  Steinborn  [19]  thermodynamically  related  the 
crystallization  pressure  to  the  supersaturation.  The  conditions  for 
reaching  equilibrium  were  obtained  by  calculating  in  two  different  ways 
the  amount  of  work  needed  to  transform  a  supersaturated  solution  into 
a  saturated  solution.  The  maximum  work  W  gained  to  perform  this  trans¬ 
formation  isothermally  and  reversibly  is 

W  =  N  kT  An  —  (2-1) 

00 

where  C  is  the  actual  concentration,  is  the  saturated  concentra¬ 
tion,  N  is  the  number  of  molecules  transformed,  k  is  the  Boltzmann's 
constant  and  T  is  the  absolute  temperature.  The  crystal  grows  under 
a  load  B  and  the  work  gained  by  lifting  the  load  for  a  distance  of  d 
is 


W  =  Bd  =  PAd  =  PN  vg  (2-2) 

where  P  is  the  crystallization  pressure,  A  is  the  contact  area  between 
the  load  and  crystal,  v$  is  the  molecular  volume  of  the  crystal. 
Equating  the  above  two  equations,  the  crystallization  pressure  is 
obtained  as 


P  = 


(2-3) 


The  above  relation  indicates  the  crystallization  pressure  in¬ 
creases  with  the  supersaturation.  Experimental  results  [4,19]  on  alum 
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covered  by  a  glass  plate  confirmed  the  relation  at  small  and  medium 
supersaturation  (C/C*,  <  1.2).  However,  the  theoretical  prediction 
was  too  large  at  high  supersaturation  [P  =  55  Kg/cm2  for  theoretical, 
while  42  Kg/cm2  was  measured  on  a  (111)  face  and  32  Kg/cm  on  a  (110) 
face].  The  (100)  face  of  alum  gave  zero  pressure  at  all  super¬ 
saturations  [19,20].  It  was  also  found  that  the  pressure  depended  on 
the  nature  of  the  external  body.  Alum  covered  by  mica  yielded  zero 
pressure  for  (111),  (HO)  and  (100)  faces  at  all  supersaturations. 
Therefore,  the  crystallization  pressure  cannot  be  described  entirely 
by  the  above  relation.  It  depends  not  only  on  supersaturation,  but 
also  on  the  type  of  crystal,  on  the  solvent,  on  the  external  body,  and 
on  the  transport  conditions  at  the  contact  area  between  the  crystal 
and  the  external  body. 

Khaimov-Mal'kov  [21]  also  studied  crystallization  pressure 
thermodynami cally .  Although  his  theoretical  relation  included  the 
anisotropy  effects  of  elastic  behavior,  these  effects  became  appre¬ 
ciable  only  at  pressures  of  geological  order  (greater  than  1000  Kg/ 

O 

cm  ). 


2.  Correns1  Phase-Boundary  Force 

As  indicated  above,  a  (111)  face  of  alum  pushes  a  glass  plate 
while  a  (100)  face  sticks  to  such  a  plate.  Further,  both  faces  stick 
to  a  mica  plate.  These  experimental  results  give  the  impression  that 
the  crystallization  pressure  might  be  related  to  surface  energies. 
Correns  [22]  established  another  condition  to  show  whether  or  not  an 
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external  body  could  be  pushed.  When  a  particle  (p)  is  in  intimate 
contact  with  a  crystal  (s)  in  liquid  (£),  an  amount  of  work  is 
necessary  to  separate  the  particle  from  the  crystal.  Work  has  to  be 
done  against  the  surface  energies  in  order  that  two  new  solid-liquid 
and  particle-liquid  interfaces  appear  and  the  original  particle-solid 
interface  disappears.  If  the  surface  energy  between  the  solid  and 
particle  is  osp,  the  surface  energy  between  the  solid  and  liquid  is 
aS£,  and  the  surface  energy  between  the  particle  and  liquid  is  ap£, 
the  work  done  is 


W  =  (a 


sp 


-  a 


S£ 


-  A 


(2-4) 


If  °sp  5  °SH  +  V  then  W  is  positive  and  work  is  gained  so  that  the 
crystal  will  lift  the  particle,  provided  that  sufficient  solution  (or 
melt)  enters  into  the  contact  region  between  the  crystal  and  the 
particle.  On  the  other  hand,  when  a  <  a  +  a  is  valid,  the 

j  P  O  X/  |JX; 

particle  simply  adheres  to  the  crystal. 

Unfortunately,  there  are  no  experimental  data  available  for 

°  .  There  are  few  values  available  for  o  „  and  a  „[23l.  The  rela- 

sp  p£  s£ 

tionship  osp  >  +  ap^  only  indicates  whether  or  not  it  is  possible 

for  the  growing  crystal  to  push  the  external  body. 


3.  Work  of  Corte 

Corte  [5]  investigated  pushing  of  seven  different  particles  by 
upward  freezing  of  water.  For  each  particle  it  was  found  that  trapping 
occurred  only  when  critical  freezing  rate  was  exceeded.  Using 


different  shapes  of  silica  (glass  bead,  broken  glass  and  quartz),  it 
was  demonstrated  that  an  important  factor  in  particle  repulsion  and 
capture  was  the  shape  of  the  particle.  Glass  beads,  having  the 
smallest  contact  area,  showed  the  smallest  freezing  rate  for  trapping 
Broken  glass  and  quartz,  having  a  greater  contact  area,  were  pushed 
at  larger  rates.  Mica  particles,  which  had  the  largest  contact  area 
of  the  particles  (all  0.15  mm  in  diameter),  exhibited  the  greatest 
critical  freezing  rate. 

The  relation  of  critical  freezing  rate  to  particle  size  was 
also  studied.  Fine  particles  migrated  under  a  wide  range  of  rates  of 
freezing,  while  coarse  particles  migrated  at  slower  and  narrower 
ranges  of  rates  of  freezing.  He  pointed  out  that  a  layer  of  water 
must  be  continuously  present  between  the  particle  and  the  ice  front 
during  particle  migration. 

4.  Work  of  Uhlmann,  Chalmers,  and  Jackson 

Uhlmann,  Chalmers,  and  Jackson  [6]  studied  the  interaction 
between  particles  and  an  advancing  solid-liquid  interface  both 
experimentally  and  theoretically.  They  observed  the  pushing  of 
eleven  different  types  of  particles  in  three  organic  compounds  con¬ 
tained  as  a  thin  film  (300-500  ym)  between  glass  slides.  They  also 
found  that  for  each  particle-melt  combination  a  critical  velocity  \1Q 
existed,  below  which  the  particles  were  rejected  by  the  interface, 
and  above  which  they  were  trapped  in  the  crystal.  The  dependence  of 
Vc  on  various  properties  of  particles  and  matrix  material  was  studied 


A  theory  was  developed  based  on  the  assumptions  of  a  particle-interface 
repulsion  (due  to  the  differed:  of  surface  energies,  AoQ  =  o$p  - 
°sl  ~  °p2^  combined  with  diffusion  of  melt  into  the  contact  region. 

The  critical  trapping  velocity  for  a  smooth  spherical  particle  on  a 
planar  interface  was  expressed  as 

Vc  =  \  (n+1 )(La  VjD/kT  R2)  (2-5) 

where  D  is  the  diffusion  coefficient,  L  the  latent  heat  of  fusion, 
v£  the  atomic  volume  of  melt,  a  the  molecular  diameter,  R  the  radius 
of  particle,  k  the  Boltzmann's  constant,  and  n  a  constant  which  is  an 
exponent  relating  variation  of  Aoq  with  separation  distance. 

They  also  treated  the  problem  in  more  detail,  including  the 
effects  of  viscous  drag  and  particle  roughness.  However,  the  results 
were  riot  satisfactory.  One  of  the  reasons  is  a  lack  of  microscopic 
properties  for  the  supercooled  thin  film  between  the  particle  and  the 
freezing  solid.  In  their  calculations,  they  had  to  assume  a  diffusion 
coefficient  400  times  smaller  than  the  value  in  the  bulk  liquid  in 
order  to  obtain  agreement  with  the  experiment  results. 

Table  1  summarizes  their  measured  Vc,  along  with  those  observed 
by  others  [5,7,9,10,24], 

5.  Work  of  Hoekstra  and  Miller 

Hoekstra  and  Miller  [7]  postulated  that  the  transport  of  water 
in  a  thin  film  between  a  glass  particle  and  ice  is  by  diffusion.  They 
took  the  free  energy  gradient  as  a  driving  force  and  related  it  to  the 
temperature  gradient,  deriving 


TABLE  1 

SUMMARY  OF  CRITICAL  VELOCITIES  FROM  PREVIOUS  STUDIES 
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2dAcp0 

v  = - L 

Vc  kT  •  pr 


(*) 


1\  dT 

37 


(2-6) 


where  d  is  the  thickness  of  the  thin  film,  A  is  the  mobility  of  the 
liquid  film,  pg  is  the  density  of  ice,  p^  is  the  density  of  water,  and 
c  is  a  constant. 

Their  experiments,  executed  like  Corte's  [5],  confirmed  the 
linear  relationship  between  Vc  and  1/R  for  both  glass  beads  and  glass 
cylinders.  The  glass  cylinders  gave  higher  values  for  Vc  than  did  the 
glass  beads.  However,  they  did  not  test  experimentally  the  dependence 


of  Vc  on  the  temperature  gradient. 


6.  Work  of  Bolling  and  Cisse 

Bolling  and  Cisse  [8,9]  combined  fluid-flow  viscous  drag  and 
diffusion  considerations  to  derive  theoretical  relations  for  Vc>  They 
also  considered  the  indentation  of  the  interface  shape  generated  by 
the  interaction  between  the  particle  and  the  interface.  Using  an 
upward  moving  interface  to  push  four  different  particles  by  freezing 
water,  their  experimental  results  [9]  supported  well  the  deduced 
relations.  They  established  a  simple  inverse  power  law  between  Vc  and 
size  for  smooth  particles.  When  the  theory  was  made  more  complex  by 
considering  the  effects  of  particle  roughness,  gravity  and  thermal 
conductivities,  the  fit  with  experiment  was  improved.  As  predicted  by 
their  theory,  the  ratio  \IQ:  /2  Vc:  /3  Vc  was  measured  for  the  critical 
velocities  when  trapping  of  the  same  type  and  size  of  particles 
Occurred  at  a  grain  surface,  a  grain  boundary  groove  and  a  grain 


boundary  triple  point. 

Their  theoretical  development  appeared  to  be  significant  for 
understanding  trapping  and  rejection  of  insoluble  particles,  but  it 
is  very  difficult  to  follow  and  understand.  Thus,  the  study  of  their 
theory  became  part  of  the  task  of  this  dissertation.  Part  of  the 
theory  relevant  to  this  work  is  explained  here  in  detail  for  a  better 
understanding  of  the  phenomena  of  particle  pushing. 

The  steady-state  interface  shape  during  the  pushing  of  a 
foreign  particle  depends  on  the  interaction  between  the  particle  and 
the  freezing  interface.  The  interface  must  be  responsive  to  all 
forces  exerted  on  the  particle  such  as  gravity  and  viscous  drag.  This 
quite  often  results  in  a  smooth  shallow  indentation  in  the  interface 
adjacent  to  the  particle.  Crystal  growth  can  be  described  by  specify¬ 
ing  the  solid-liquid  interface  temperature  as  a  function  of  position. 
For  the  simplest  case  of  pure  growing  material  and  an  inert  smooth 
particle  of  the  same  thermal  conductivity  as  the  solid  and  melt,  the 
interfacial  temperature  Ti  can  be  expressed  as 

T,  -  -  ATs  -  ATc  -  4Ti  (2‘7 

where  T  is  the  melting  temperature  of  pure  material,  ATg  is  the 
supercooling  needed  to  provide  a  driving  force  for  crystal  growth  and 
is  a  function  of  growth  rate  V,  ATc  =  oK/AS  is  the  temperature  depres¬ 
sion  due  to  curvature  effects,  AS  is  the  entropy  of  fusion  per  unit 
volume,  K  is  the  total  curvature  at  any  point  on  the  interface,  and 
ATi  is  the  temperature  change  due  to  interaction  of  the  particle  with 
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the  interface  and  is  evaluated  as  follows. 

The  thermodynamic  criteria  for  equilibrium  requires  that  the 
free  energies  of  the  two  phases  (solid  and  liquid)  are  equal  at  the 
melting  point.  The  fundamental  relation  for  a  small  change  in  free 
energy  of  the  system  dG  is  expressed  as 


dG  =  -SdT  +  V  dP  (2-8) 

where  S  is  the  entropy  of  system,  V  the  volume,  and  P  the  pressure. 

It  follows  that  the  condition  for  equilibrium  between  two  phases  is 
that  dG$  =  dG^,  or 

-Ss  dTs  ♦  Vs  dPs  =  -St  dT,  *  »t  dPt  (2-9) 

where  the  subscripts  s  and  i  refer  to  solid  and  liquid  respectively. 

If  we  assume  that  the  liquid  pressure  in  the  particle-interface  con¬ 
tact  area  is  similar  to  the  bulk  pressure,  then  the  dP^  may  be 
neglected.  Further,  Ts  =  T^  at  the  equilibrium  state.  Thus,  the 
above  equation  becomes 


— -V-  AP  .  (2-10) 

s  ij 

If  we  substitute  AS  =  Ss~ S^/ Vs ,  and  AP  =  F/fir^  into  the  equation,  the 
temperature  change  due  to  a  force  acting  on  the  crystal  is 


4Ti  - 


ASirr, 


if  r  <  r 
—  o 


if  r  >  rQ 


(2-11) 
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where  rQ  is  an  effective  contact  radius  beyond  which  no  interaction 
occurs  between  the  particle  and  the  interface,  as  defined  by 
Uhlmar.n  et  al .  [6].  F  is  the  interacting  force  within  r  . 

The  search  for  a  complete  solution  of  F  is  the  next  step  in 
development  of  the  theory.  The  viscous  drag  on  the  particle,  which 
produces  a  force  on  the  freezing  interface,  is  considered  next.  A 
hydrodynamic  approach  is  used  to  solve  the  problem  of  flow  around  a 
sphere  to  a  nearby  uniform  sink  (interface).  Figure  1  shows  a 
schematic  cross  section  of  a  spherical  particle  in  contact  with  an 
interface.  The  contact  region  within  rQ  is  expanded.  The  formula 
for  viscous  drag  Fc  for  the  plane  sink  (0  =  0)  as  derived  in  Uhlmann's 
dissertation  by  G.  Carrier  [25]  is 

Fc  =  6-rrnV  R2/d  (2-12) 

which  gives  the  viscous  drag  upon  a  sphere  of  radius  R  located  a  dis¬ 
tance  d  from  a  plane  interface  advancing  at  growth  rate  V.  The 
problem  considered  here  requires  a  sink  which  is  concave  due  to  inter¬ 
action  within  the  contact  area  r  <_  r  .  The  separation  distance  h 
between  particle  and  interface  at  any  point  in  the  contact  region  is 

h  =  d  +  R  -  (R2-r2)1/2  -  E(r)  (2-13) 

where  E(r)  is  the  vertical  distance  between  a  concave  interface  and 
its  lowest  point,  and  is  a  function  of  the  radial  distance.  For  the 
case  of  r  <<  R,  (R  - r  )  'is  approximately  equal  to  R  -  r  / 2R  by  a 
MacLaurin's  series  and  thus 
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Figure  1  Schematic  cross  section  of  a  spherical 
particle  in  contact  with  a  solid-liquid 
interface  with  a  magnified  contact  area. 
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h  »  d  +  r2/2R  -  E(r)  .  (2-14) 

A  parabolic  velocity  profile  U  is  assumed  for  the  flow  in  the  contact 
region,  with  zero  tangential  velocity  at  the  solid-liquid  interface 
(z  =  E(r))  and  the  particle-liquid  interface  (z  =  h  +  E(r)).  That  is, 

U  =  C(r)  (z  -  E(r)j  (z  -  h  -  E(r ))  .  (2-15) 

Consider  the  region  of  contact  as  a  sort  of  cylinder  with  a 
radius  of  r.  Conservation  of  mass  requires  that  the  liquid  entering 
the  cylindrical  walls  equals  that  leaving  at  the  bottom  of  the 
cylinder  under  the  steady  state  conditions. 

r  h+E(r )  2 

2tt  rU  p  dz  =  nr  V  p  .  (2-16) 

'  E(r)  *  S 

Substituting  Eq.  (2-15)  into  the  above  equation  and  integrating  the 
left-hand  side,  the  function  C(r)  is  found  to  be 


Hence,  the  complete  expression  for  U  is 

u  =  f^fe)(z'E(r))(z'h'E(r)^  '  l2",8) 

The  pressure  gradient  in  the  r  direction  is  found  from  the 
Navier-Stokes  equation  for  a  fluid  in  steady-state  laminar  motion  [26] 


to  be 
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iE  =  n  A  =6arV/V) 

dr  V  h3  W  '  ‘ 

The  pressure  difference  AP  between  the  inside  of  the  contact  region 
and  the  outside  of  the  contact  region  is 

r 

R 

Since  h  increases  with  r  in  real  situations,  dP/dr  in  Eq.  (2-19) 

increases  from  zero  and  then  decreases,  approaching  zero  again  as  r 

approaches  R.  On  the  other  hand,  AP  is  the  largest  at  r  =  0  and 

decreases  as  r  increases.  Their  variations  are  shown  schematically 

in  Figure  2.  These  trends,  as  given  by  Eqs.  (2-19)  and  (2-20),  hold 

whether  the  approximate  formula  Eq.  (2-14)  is  employed  or  not.  In 
3 

either  case  r/h  approaches  zero  as  r  approaches  R.  The  major  con¬ 
tribution  to  the  integral  Eq.  (2-20)  comes  from  the  region  r  <<  R. 

Because  the  hydrostatic  pressure  P(R)  can  be  considered  con¬ 
stant  over  the  surface  of  a  small  particle,  AP  is  essentially  the  net 
pressure  exerted  on  the  particle  in  the  vertical  direction.  This 
pressure  is  multiplied  by  the  projected  area  of  the  particle,  and 
thus  gives  a  net  force  F,  which  is 


Since  we  have  assumed  r  «  R,  we  may  well  take  the  lower  limit  as 


AP  =  P(r)  -  P(R)  = 


(2-21) 


Particle 


Figure  2 


£P 


Sketch  showing  variations  of  dP/dr  and  ^P 
with  respect  to  r  in  the  contact  region^ 
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00  for  R  to  obtain 


The  above  equation  cannot  be  integrated  unless  we  know  the  function 
E(r).  Because  we  are  searching  for  an  outer  bound,  it  is  possible  to 
find  a  limit  for  this  function  which  may  give  a  maximum  force.  At 
the  limit  r  =  0  and  z  =  0,  the  contact  region  between  the  particle  and 
interface  may  approximate  a  part  of  a  sphere,  and  E(r)  has  a  maximum 
value  if  the  interface  would  closely  contact  with  a  spherical 
particle  with  a  radius  of  R/a  where  a  £  1.  At  the  limit  d  «  h,  we 
obtain  from  Eq.  (2-14) 


E(r) 


a  r 
~2R 


(2-23) 


Equation  (2-22)  is  then  integrated  by  substituting  the  above 
relation  for  C(r)  to  obtain 


F(g) 


S5  6Trn  V  R2  / £s\ 
d(l-g)2 


(2-24) 


where  F(a)  £  F. 

When  a  has  a  value  of  zero,  the  interface  is  a  plane  sink  and 
not  responsive  to  the  presence  of  the  particle.  Equation  (2-24)  thus 
becomes  Eq.  (2-12)  of  Carrier's  expression  for  viscous  drag.  At  the 
other  extreme,  when  a  has  a  value  of  unity,  the  interface  and  the 
particle  are  in  an  intimate  contact,  and  Eq.  (2-24)  yields  an 
extremely  large  force.  Therefore,  we  have  to  limit  the  use  of  F(g), 
particularly  as  a  approaches  unity. 
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The  viscous  drag  derived  above  appears  to  depend  on  the 
separation  distance  d  at  r  =  0.  The  value  of  d  is  difficult  to 
evaluate  so  that  we  may  relate  it  to  some  system  parameters.  The 
interface  discussed  here  is  responsive  to  the  existence  of  the 
nearby  particle;  and  therefore  d  has  an  influence  on  the  interface 
shape.  If  we  relate  the  interface  shape  parameter  a  to  some 
characteristic  value  of  d,  let  us  say  at  the  value  of  NaQ  where  aQ 
is  the  interatomic  distance  and  N  is  a  number  depending  on  a  system 
parameter  3  (3  is  some  characteristic  only  of  the  liquid  for  inert 
particle  but  of  liquid  and  particle  if  there  is  a  chemical  inter¬ 
action),  and  choose  a  decay  function  to  represent  the  variation 
between  the  two  limits  for  a  (a  =  0  and  1),  we  obtain 


a  =  exp  [- (d-3aQ )/ aQ]  ,  (2-25) 

which  gives  a  planar  interface  (a  =  0)  at  d  »  3aQ  and  intimate  con¬ 
tact  (ct  =  1)  between  the  particle  and  the  interface  at  d  =  3aQ. 
Rearranging  the  above  relation  we  obtain 


d  =  aQ{3-i>n  a)  .  (2-26) 

Thus,  Eq.  (2-24)  can  be  expressed  in  terms  of  3  by  means  of  Eq.  (2-26) 


6tt  n  V  R2  /M 
(l-cx)2(3-£n  a)aQ  \pfc' 


(2-27) 


We  just  completed  the  hydrodynamic  approach  to  the  transport 
problem  in  the  contact  region.  Uhlmann  et  al.  [6]  treated  the  same 
transport  problem  from  a  diffusion  viewpoint.  If  there  is  a  diffusion 
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limit,  Eq.  (2-24)  may  be  used  to  compare  with  this  diffusion  limit. 

For  the  purpose  of  comparison,  F(a)  is  mathematically  manipulated  with 

o 

the  assumption  of  w  =  r  ( 1  -ot)/2R , 


F(a)  =  - 


6tt  n  V  R2  /  ps 

\  f°  dw 

(1-ct)2  W' 

/  -O  (d+w)2 

6tt  n  V  R2  /V 

)\\ 

fr°  dw 

(1  -a)2  W' 

'u 

00  (d+w)( 

dw 


rQ  (d+w) 


(2-28) 


The  force  F(a)  is  divided  into  two  portions.  One  arises  within 
the  contact  region  r  £  r  ,  the  other  one  from  the  region  r  >  rQ.  Now 
we  may  compare  F(a)  with  a  supposed  force  which  is  required  for  the 
diffusion  process , 


virtual 


_  _  6tt  n  V  R 


(1-ct) 


2 


•G)J 


dw 


(d+w) 


2  +  FD 


(2-29) 


where  Fp  is  the  force  needed  for  diffusion  in  the  contact  area 
r  £  rQ.  This  force  is  less  than  that  for  flow  within  the  same 

region,  or  flow  would  occur.  Thus,  F(a)  £  ^virtual’  ^  we  wor*<  at 
the  limit,  it  becomes  immaterial  to  distinguish  between  flow  and 
diffusion.  The  movement  of  molecules  on  the  surface  following  a  ran¬ 
dom  walk  with  no  preferred  direction  is  [27] 

Ar2  =  4  Dt  (2-30) 


where  t  is  the  time  necessary  for  a  mean  displacement  Ar.  If  we 
assume  that  b  is  the  height  for  growth  of  a  monoatomic  layer,  then 
the  time  required  for  molecules  to  di  rfuse  into  the  contact  region 


mM 


(with  an  effective  radius  r  )  is  t  =  b/V.  The  least  advantageous 
diffusion  sequence  is  then 


<  4  D  b/V 


(2-31) 


If  we  relate  fluid  flow  and  diffusion  by  means  of  the  Stokes-Einstein 


equation,  that  is,  D  =  kT/37rnao,  the  above  relation  becomes 

r2  <  JJJ„ 

ro  -  3tt  n  V  • 


(2-32) 


With  this  relation  in  mind,  we  may  determine  the  conditions 
under  which  the  force  F(a)  is  dominant  during  particle  pushing.  The 
force  within  r  <  rQ  where  the  diffusion  might  occur  can  be  evaluated 
from  the  second  term  of  Eq.  (2-28) 


'o  to  r. 


in  r)  V  /Ps  \  f° 
(1-a)^  '  ■*  r 

n  V  R2  /Ps\ _ 

1-rvl  \p!./  ?RHf 


rQ  (d+w)‘ 


6m  n  V  R  ro _ 

(1  -a)  W'  2Rd[d+(l-a)r‘"/2R] 


(2-33) 


The  force  F(a)Q  tQ  r  is  sensibly  the  force  F(a)  of  Eq.  (2-24)  if 

o 


rg(l-a) 


>>  d  . 


(2-35) 


At  the  limit  of  the  diffusion  sequence,  the  above  relation  when  com¬ 
bined  with  Eq.  (2-32)  becomes 


3tt  n  V  R 


>>  d  . 


(2-35) 
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Based  on  Eq.  (2-24),  the  viscous  drag  increases  with  increase 
of  the  growth  rate  or  the  interface  curvature.  F(a)  reaches  a 
maximum  when  the  growth  is  at  the  critical  trapping  velocity.  At 
this  point  the  region  of  close  contact  gives  a  maximum  radius  R/a. 

If  we  equate  the  effect  of  the  maximum  force  to  that  of  the  maximum 
possible  curvature,  or  ATc  =  AT^,  we  obtain  the  conditions  for 
initiation  of  particle  trapping.  Thus 


F(a)  2o  a 
AC  T  =  AS  R 
AS  "  ro 


or 


3n_y; 

aa_ 


1 _ 

a(l-a)  (&-£n  a) 


1  • 


(2-36) 


If  Eq.  (2-35)  holds,  we  are  allowed  to  use  Eq.  (2-32)  to 
7 

evaluate  the  quantity  rQ  while  maintaining  the  conditions  of  maximum 

force  and  critical  trapping  velocity.  That  is, 

4  kT  oa. 

1  ♦<«>  ( i 

s 


- * 

9tt  n  R 


(2-37) 


where  <t>(ot)  is  a(l  ot)2(B-£n  a).  This  is  the  final  equation  to  be  used 
for  prediction  of  the  critical  trapping  velocity. 

The  relation  of  a  and  8  is  found  at  the  extremum  of  <|>(a),  or  at 
the  condition  of  d4>(a)/da  =  0.  They  are 


6  =  £n  a  +  (l-a)/(l-3a) 


(2-38) 


tWBI  *  a(l-a)3/(l-3a) 


and 


(2-39) 
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The  above  relations  restrict  the  quantity  a,  which  must  be 
less  than  1/3  or  else  the  curvature  of  the  interface  would  be  large 
enough  to  engulf  the  particle.  If  we  apply  this  limitation  to  Eq. 
(2-25),  a  useful  relationship  is  obtained  as 


or  ■ 

d  >  (1.1  +  3)aQ  .  (2-40) 

i  } 

The  above  relation  gives  the  thickness  of  the  thin  film 
existing  between  the  particles  and  interface.  Since  the  thin  film 
consists  of  several  interatomic  layers,  it  is  logical  to  assume  a 
minimum  of  unity  for  6  to  solve  the  problem. 

The  relationship  between  d  and  a  can  be  obtained  from  Eqs. 

i 

(2-26)  and  (2-38), 

I 

d  =  (T-a)a0/(l-3a)  .  \  (2-41) 

,  |  I 

The  inequality  presented  in  Eq.  (2-35)  to  be  evaluated  then 

i  , 

becomes  , 

VR  «  2(l-3a)kT/ 3tt  n  aQ  .  (2-42) 

i  '  i 

If  we  choose  6  =  1 ,  the  value  of  a  gives  0.229  at  4>(ot)  = 

-8 

iiia)  .  For  water  as  a  matrix  material ,  a  =  3.8  x  10  cm, 
n  =  1.8  c. p. ,  T  =  273°K,  and  VR  <  3.7  x  10'6  cm2/sec:  This  require¬ 
ment  was  met  by  several  different  particles  in  freezing  of  water  [7,9]. 
Equation  (2-27)  was  derived  without  accounting  for  effects  of 

•i  i  , 

gravity  and  particle  roughness,  Cisse  and  Bolling  [9]  studied  smooth 


I 


! 


I 
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spherical  tungsten  and  copper  particles  in  water.  The  data  supported 
Eq.  (2-37)  quite  well  for  the  small  particles  with  R  <  10  pm. 

B.  Some  Theoretical  Considerations 

1 .  Structure  of  the  Solid-Liquid  Interface 

The  importance  of  interface  morphology  to  this  study  is 
obvious,  since  the  solid-liquid  interface  was  employed  to  reject  and 
push  foreign  particles.  The  nature  of  the  interface  is  known  to  have 
a  decisive  influence  on  the  kinetics  and  morphology  cf  crystal  growth. 
According  to  one  viewpoint,  all  molecules  either  clearly  belong  to  the 
crystal  or  to  the  melt,  yielding  a  sharp  well-defined  interface 
between  the  two.  Such  an  interface  may  either  be  atomically  smooth 
or  atomically  rough.  These  two  possible  interfaces  are  shown  in 
Figure  3  where  the  white  spheres  represent  the  atoms  of  the  melt,  the 
black  spheres  are  the  crystal  atoms.  Figure  3(A)  is  an  atomically 
smooth  interface,  while  (B)  represents  a  rough  interface. 

With  the  present  knowledge  of  crystal  growth,  two  possible 
growth  mechanisms  have  been  discussed  [28,29].  A  crystal  whose 

1 

equilibrium  interface  is  atomically  smooth  grows  by  the  formation 
and  lateral  growth  of  layers.  In  a  perfect  crystal,  new  layers  are 
formed  discontinuously  by  two-dimensional  nucleation.  When  one  or 
more  screw  dislocations  are  present,  new  layers  may  be  formed  con¬ 
tinuously.  Rapid  lateral  growth  of  layers  leads  to  formation  of 
macroscopic  facets. 

1  When,  on  the  other  hand,  the  equilibrium  interface  is  atomically 


(A)  Atomically  smooth  interface 


Crystal 


(B)  Atomically  rough  interface 


Crystal 


Figure  3  Atomic  solid-liquid  interface 
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rough,  the  crystal  can  grow  without  formation  of  new  layers.  The 
interface  advances  uniformly  and  is  not  faceted. 

The  theory  of  the  equilibrium  interface  roughness  was  satis¬ 
factorily  examined  by  Jackson  in  1958  [30].  He  considered  an  initially 
atomic  smooth  interface  and  calculated  the  change  in  surface  free 
energy  on  adding  extra  molecules.  The  results  of  his  elementary 
statistical  thermodynamic  approach  predicted  that  the  equilibrium 
interface  should  be  atomically  rough  and  non-faceted  for  materials 
having  entropies  of  fusion  (aS^)  less  than  2R  or  16.8  J/mole°K  (where 
R  is  the  gas  constant).  For  ASf  >  4R  the  interface  was  predicted  to 
be  atomically  smooth,  faceted  and  essentially  crystal lographical ly 
perfect.  Most  metals  and  certain  organic  compounds  have  low  entropies 
of  fusion  (ASf  <  16.8  J/mole°K).  Their  non-faceted  and  planar  inter¬ 
faces  were  confirmed  by  recent  experiments  [29,31]. 

A  different  approach  to  the  problem  of  the  structure  of  the 
interface  was  discussed  by  Cahn  [32]  who  assumed  a  diffuse  interface 
of  finite  thickness.  The  transition  from  crystal  to  melt  takes  place 
over  several  layers  of  atoms.  According  to  the  theory,  lateral 
growth  occurs  at  low  undercooling,  but  continuous  growth  takes  place 
at  high  undercooling.  A  transition  from  lateral  growth  kinetics  to 
continuous  growth  kinetics  was  predicted  for  all  materials.  Unfor¬ 
tunately,  there  is  no  known  way  of  determining  the  diffuseness  of  the 
interface,  so  that  a  test  of  the  theory  is  impossible. 

2.  Interfacial  Temperature  Gradient 

The  temperature  gradient  at  the  sol  id- liquid  interface  is  an 
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important  parameter  in  the  growth  of  crystals.  It  influences  the 
perfection  and  propertie..  of  crystals  [33],  eutectic  structures  [34], 
the  microscopic  shape  stability  of  the  interface  [35],  and  interface 
faceting  [36].  The  role  of  temperature  gradient  in  pushing  and 
trapping  is  not  clear.  Hoekstra  and  Miller  [7]  predicted  that  Vc 
should  be  proportional  to  the  temperature  gradient.  However, 
theoretical  treatments  by  other  workers  [6,8]  did  not  consider  this 
term.  The  importance  to  this  work  is  discussed  below. 

Interface  Stability.  Even  so-called  "pure"  materials  are  not 
completely  pure.  Traces  of  impurities  are  always  present.  If,  for 
example,  the  concentration  of  the  solute  (impurity)  in  the  solid  is 
less  than  that  of  the  liquid  from  which  it  is  forming,  there  is  a 
rejection  of  solute  into  the  liquid  at  the  solid-liquid  interface. 

The  liquid  in  the  vicinity  of  the  interface  becomes  enriched  in 
solute.  Unless  the  growth  rate  is  slower  than  the  diffusion  rate  of 
solute  in  the  liquid,  a  concentration  gradient  develops  in  the  liquid 
ahead  of  the  interface.  Since  the  presence  of  solute  in  the  melt 
lowers  the  freezing  points,  the  freezing  temperature  increases  from 
the  interface  while  the  impurity  concentration  decreases  as  shown  in 
Figure  4.  The  imposed  temperature  normally  increases  as  one  moves 
into  the  melt.  When  the  slope  of  this  actual  temperature  is  less  than 
the  slope  of  the  freezing  point,  the  melt  ahead  of  the  interface  is 
below  its  normal  freezing  temperature.  This  excess  supercooling, 
termed  constitutional  supercooling  [37],  is  a  direct  result  of  the 
concentration  gradient  which  exists  in  the  melt.  Under  such  conditions 


Concentration 
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(A)  Solute  distribution  during  normal  freezing. 


■^Distance  into  melt 


(B)  Corresponding  equilibrium  freezing  temperature 
distribution. 
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Figure  4  Constitutional  supercooling. 
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the  interface  is  unstable  with  respect  to  small  perturbation.  The 
supercooling  can  be  relieved  by  the  growth  of  small  projections  into 
the  melt.  A  freezing  interface  which  is  rough  on  an  atomic  scale  may 
develop  pox  structures,  elongated  cells,  hexagonal  cellular  structures, 
or  even  dendrites,  depending  on  the  degree  of  constitutional  super¬ 
cooling.  Such  interface  breakdown  has  a  detrimental  influence  on  the 
process  of  particle  pushing  and  therefore  was  avoided  in  this  work. 

The  conditions  for  the  onset  of  constitutional  supercooling  may 
be  expressed  in  simple  mathematical  terms  [38,39].  When  the  imposed 
temperature  gradient  is  less  than  the  gradient  of  equilibrium  freezing 
temperature  at  the  interface,  the  following  relationship  between  the 
imposed  growth  conditions  and  the  system  parameters  is  obtained 
for  the  critical  condition  of  supercooling  if  there  is  no  liquid 
mixing  in  the  liquid  phase. 


<  m 


<cs>i  -B7(ct>i 


(2-43) 


where  (C  ).  and  (Cj-  are  the  interfacial  solute  concentrations  in 
the  solid  and  the  melt,  and  m  is  the  slope  of  the  liquidus  line  on 
the  phase  diagram.  Some  quantitative  studies  [40,41]  have  verified  trie 
above  relationship.  In  order  to  avoid  the  effects  of  constitutional 
supercooling,  the  imposed  temperature  gradient  should  be  increased  to 
exceed  the  gradient  of  freezing  temperature  at  the  interface. 


Interface  Faceting.  Jackson's  theory  of  interface  roughness 
was  derived  for  an  isothermal  and  equilibrium  environment.  In  real 
solidification  situations,  the  system  is  neither  isothermal  nor  at 
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equilibrium.  Wilcox's  theoretical  discussions  [36]  pointed  out  that 
crystals  predicted  to  grow  in  a  faceted  manner  (high  entropy  of  fusion 
materials)  would  grow  non-faceted  in  the  presence  of  a  sufficiently 
large  temperature  gradient.  From  equilibrium  cons’H' rations ,  faceting 
in  a  temperature  gradient  would  exist  only  if  the  edge  energy  was 
considered  in  addition  to  the  surface  energy  and  the  bulk  free  energy. 
From  the  standpoint  of  kinetic  effects,  faceting  is  favored  by 
anisotropic  growth  kinetics  and  rapid  growth.  There  ir  widespread 
qualitative  experimental  evidence  to  support  the  effect  of  temperature 
gradient  on  interface  faceting. 


CHAPTER  III 


EXPERIMENTAL 

Two  crystal  growth  methods  were  employed  for  this  experimental 
work.  A  vertical  Bridgman-Stockbarger  technique  was  first  used  to 
investigate  the  fundamental  nature  of  particle  pushing  and  trapping 
phenomena.  Horizontal  zone- ref i ning  with  rotation  was  then  em¬ 
ployed  to  study  some  of  the  design  variables  for  a  separation  process. 
Uhlmann,  et  al .  [6]  conducted  particle  pushing  experiments  by  means  of 
a  horizontal  thin  film  (held  between  glass  slides)  with  a  vertical 
freezing  interface.  Their  results  might  have  been  strongly  influ¬ 
enced  by  the  glass  slides  and  not  characteristic  of  bulk  materials. 
Others  [5,7,9,10]  used  Bridgman-Stockbarger  growth  with  the  solid- 
liquid  interface  moving  upward,  reducing  the  possibility  of  the 
particles  dragging  on  the  glass  wall.  We  incorporated  facilities  for 
direct  observation  and  temperature  measurements.  This  was  the  first 
time  that  a  horizontal  zone-refiner  with  rotation  was  used  to  investi¬ 
gate  removal  of  particles  from  the  matrix  materials. 

Following  the  description  of  the  two  experimental  techniques, 
the  materials  used  for  this  study  are  presented,  including  the 
particles  and  the  organic  compounds.  The  treatment  of  solidified 
products  is  then  discussed.  The  final  section  of  this  chapter  is 
devoted  to  a  description  of  the  zone-refining  of  organic  compounds. 
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A.  Bridgman-Stockbarger  Technique 
1 .  Apparatus 

A  Bridgman-Stockbarger  apparatus  was  constructed  as  shown  in 
Figure  5.  This  unit  was  designed  to  provide  an  upward  moving  solid- 
liquid  interface  which  could  be  constantly  observed  with  a  micro¬ 
scope.  Essentially  the  unit  consisted  of  a  heater,  a  cooler  and  a 
tube  lowering  mechanism.  The  heater  was  made  by  winding  ni chrome 
wire  on  a  piece  of  flanged  Pyrex  glass  tube,  which  was  insulated  by 
another  glass  tube.  The  cooler  was  constructed  similar  to  a 
laboratory  condenser,  as  shown  in  Figure  6.  The  cooling  fluid  was 
circulated  from  a  refrigerated  constant  temperature  bath  (0  to  -30°C), 
passed  through  the  annular  space  of  the  cooler,  and  returned  to  the 
circulator.  A  crystal-growth  tube  travelling  through  the  heater 
steadily  displaced  the  heat-transfer  fluid  which  overflowed  from  a 
side  tube.  This  fluid  was  placed  in  the  inner  tube  of  the  cooler  to 
eliminate  the  air  gap  between  the  inner  tube  and  the  crystal -growth 
tube,  thus  improving  the  heat  transfer.  The  circulating  fluid  and 
the  heat-transfer  fluid  were  not  mixed,  but  both  consisted  of  50 
volume  %  ethylene  glycol-water.  The  position  of  the  side  tube  was 
adjustable,  so  that  the  level  of  the  fluid  in  the  inner  tube  could  be 
controlled.  Two  variable-speed  motors  provided  lowering  rates  of 
1  to  15  mm/hour  and  10  to  80  mm/hour.  The  heater  power  was  controlled 
by  an  auto-transformer  which  was  connected  to  a  Sola  constant- 
voltage  transformer  (30  VA,  Harmonic  Neutralized  Type  CVS)  to  obtain 
a  constant  voltage  A.C.  input.  The  power  (3  to  12  watts)  was 


Figure  5  Bench-scale  unit  of  improved  Bridgman- 
Stockbarger  technique. 

(1)  Drive  mechanism 

(2)  Heater  and  cooler  (Details  shown  in  Fig. 6) 

(3)  Constant  temperature  circulator 

(4)  Wattmeter 

(5)  Electronic  controlled  stirrer 

(6)  Microscope 

(7)  Movie  camera  equipped  with  time-lapse 
controller 


Crystal  growth  tube 
traveling  through  the 
heater  and  cooler 


Figure 
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Circulating  cooling  fluid 


6  Schematic  diagram  of  the  heater 
and  the  cooler. 


monitored  by  a  wattmeter  placed  between  the  heater  and  the  auto¬ 
transformer. 
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Most  runs  were  performed  with  a  cooler  made  of  brass  (37  mm 
I.D.  x  17  cm  long  for  the  inner  tube)  while  a  few  runs  were  made  with 
a  small  glass  cooler  (17  mm  I.D.  x  15  cm  long  for  the  inner  tube). 

Two  heaters  were  employed--one  small  (13.5  mm  I.D.  x  12  cm  long  for 
the  inner  tube)  and  one  large  (31  mm  I.D.  x  14  cm  long  for  the  inner 
tube).  The  large  heater  was  used  for  some  runs  with  a  glass  stirrer 
inserted  into  the  growth  tube.  The  stirring  motor  was  placed  on  a 
adjustable  stepless  platform  in  order  that  the  distance  between  the 
stirring  blade  and  the  freezing  interface  could  be  controlled. 

2.  Sample  Preparation 

The  solidification  tube  was  made  from  Pyrex  glass  tubing  with 
a  weight  attached  to  the  round  bottom  to  keep  the  growth  tube 
vertical.  Most  runs  were  made  with  10.5  mm  I.D.  (13  mm  O.D.)  x 
10  and  15  cm  long  tubes.  Some  runs  with  stirring  used  large  tubes, 
26.5  mm  I.D.  (30  mm  O.D.)  x  22  cm  long.  In  the  early  experiments, 
particles  were  directly  added  into  an  organic  melt  and  the  mixture 
then  poured  into  the  growth  tube.  As  described  in  Section  C,  these 
particles  were  agglomerates.  They  should  have  been  broken  into 
individual  particles  before  using  for  particle  pushing  study.  In 
later  experiments  the  particle-melt  mixture  was  first  poured  into  a 
10  cc  cylindrical  Pyrex  bottle,  sealed  with  a  screw  cap,  and  then 
placed  in  an  ultrasonic  cleaner  containing  hot  water  in  order  to  keep 
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the  organic  molten.  The  mixture  was  ul trasonically  vibrated  about  2 
minutes,  and  then  transferred  into  the  growth  tube  by  a  dropping  pi  pet. 

Molten  organic  chemicals  dissolve  gases  and  release  them  as 
bubbles  at  the  freezing  interface  during  solidification.  Th^se  bubbles 
were  found  to  seriously  interfere  with  the  particle  pushing  process. 
Bubbles  were  avoided  in  later  experiments  by  the  following  procedure: 

A  pressure  slightly  less  than  one  atmosphere  was  applied  to  the 
organic  melt  in  a  tube  so  as  to  gently  degas  the  melt.  The  organic  was 
then  solidified  and  a  vacuum  (50  Torr)  applied  to  the  tube,  which  was 
sealed  with  a  torch.  Gas  bubbles  were  not  generated  during  solidifica¬ 
tion  in  these  vacuum-sealed  tubes.  Gas  bubbles  are  further  discussed 
in  Chapters  IV  and  VI. 

3.  Experimental  Procedure 

A  solidification  tube  containing  a  mixture  of  particles  and 
organic  compound  was  placed  in  the  growth  unit  as  illustrated  in 
Figure  6.  The  length  of  the  tube  submerged  in  the  cooler  was 
initially  about  one  tube  diameter.  The  reason  for  this  is  described 
in  the  next  chapter.  As  soon  as  a  steady  solid-liquid  interface 
covered  with  a  layer  of  particles  was  established,  the  tube  was  lowered 
at  a  predetermined  rate.  The  interface  was  observed  through  a  stereo¬ 
scopic  microscope  at  10-50  X  and  the  interface  position  determined 
by  the  eyepiece.  Each  run  started  with  a  low  lowering  rate  (about 
6  mm/hour)  which  was  increased  step  by  step  (about  2  mm/hour  every 
30  minutes)  until  all  particles  were  trapped  by  the  growing  solid. 
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For  some  time  after  each  increase  in  lowering  rate  the  growth  rate  was 

less  than  the  lowering  rate.  Relative  to  the  observer,  the  interface 

!  1 

moved  slowly  downward  after  a  rate  change  until  a  new  steady-state 

'  :  *  I 

position  was  reached  (stationary  with  respect  to  the  heater  and 
cooler).  Under  this  steady-state  condition,  the  growth  rate  and  the 

j 

lowering  rate  were  equal.  Normally,  this  transient  period  was  20  to 
30  minutes  in  length.  The  minimum  growth  rate  for  particle  trapping, 
V  was  usually  determined  by  direct  observation  of  pushing  and 
trapping,  but  sometimes  it  was  estimated  by  microscopic  examination 

of  the  solidified  product.  The  errors  in  such  measurements  were  less 
than  10%.  i  i  i 

Motion  pictures  were  taken  of  interesting  phenomena  at  the 
freezing  interface,  such  as  bubble  formation  and  surface-driven  flows 

I  I  I  I  •  s 

using  a  16  mm  movie  camera  on  the  microscope,  as  shown  in  Figure  5. 
Time-lapse  cinematography  was  employed  for  very  slow  particle  pushing 
processes.  To  avoid  large  disturbance  of  the  heat  balances  at  the 

■  i  ,i 

advancing  interface,  a  water  cell  was  placed  in  front  of  the  500-watt 
photo  lamp.  Water  flowed  through  the  cell  continuously  and  absorbed 
heat  generated  from  the  lamp.  The  time-lapse  movie  showed  that  the 
travelling  motion  of  the  growth  tube  was  not  smooth,  presumably 
because  of  backlash  in  the  gearing.  Additional  weight  was  placed 

i 

between  the  motor  and  the  tube  to  reduce  this  problem. 

4 •  Temperature  Measurements 

Thin  36-gauge  copper  and  constantan  thermocouple  wires  were 
employed  to  measure  the  temperature  distribution  in  the  Bridgman 
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growth  of  naphthalene.  These  thin  wires  were  supported  by  e  capillary 
tube  and  positioned  at  the  center  of  the  cross  section  of  the  growth 
container  as  shown  in  Figure  7(A).  The  thermocouple  emf  was  measured 
by  a  millivolt  recorder  against  an  ice  bath  junction. 

Because  the  thermocouple  was  fixed  in  the  growth  container, 
it  travelled  with  the  container  as  the  latter  passed  through  the 
growth  unit  (Figure  S).  Therefore,  temperature  measurements  started 
in  the  melt,  reached  the  sol  id- liquid  interface,  and  finished  in  the 
solid.  The  temperature  versus  time  was  continuously  recorded.  The 
rate  of  change  of  temperature  divided  by  the  growth  rate  yielded  the 
temperature  gradient.  Thus,  for  measuring  temperature  gradients  at 
the  interface,  knowledge  of  the  time  at  which  the  thermocouple  was 
at  the  interface  was  required.  For  instance,  if  the  thermocouple 
joint  was  0.5  mm  in  diameter  and  the  growth  rate  was  10  mm/hour,  it 
took  3  minutes  for  the  interface  to  pass  by  the  thermocouple.  This 
made  it  difficult  to  determine  the  interface  from  the  temperature 
recording,  particularly  for  experiments  at  very  slow  growth  rates. 

With  the  help  of  microscopic  observation,  the  interface  position  was 
determined  when  the  thermocouple  was  contacting  the  interface.  The 
errors  in  the  temperature  measurements  were  less  than  0.1°C,  but  they 
might  be  up  to  20%  in  the  temperature  gradients. 

Another  factor  influencing  the  accuracy  of  the  above  measure¬ 
ments  was  heat  conduction  down  the  thermocouple  wires,  so  that  the 
measured  temperatures  were  higher  than  they  should  have  been.  An 
improved  measuring  technique  was  therefore  developed  as  shown  in 


Figure  7(B).  The  bare  thermocouple  wires  were  stretched  across  the 
cross  section  of  the  container  with  the  welded  junction  positioned 
at  the  center.  The  wires  were  passed  through  two  opposing  small 
holes  in  the  container  tube  and  the  holes  were  sealed  with  epoxy 
cement.  This  improved  method  was  not  actually  used  in  the  present 
study  because  preliminary  results  based  on  the  method  of  Figure  7(A) 
indicated  that  the  critical  velocity  Vc  does  not  depend  on  the 
temperature  gradient.  However,  the  method  was  used  in  this  laboratory 
by  Chong  E.  Chang.  The  results  were  quite  satisfactory  even  at  very 
slow  growth  rates,  less  than  1  mm/hour. 

B.  Horizontal  Zone  Refining  with  Rotation 

Some  experiments  were  made  in  the  vertical  Bridgman  apparatus 
by  inserting  a  stirrer  in  the  growth  tube  (Figure  6).  The  purpose  of 
stirring  of  the  melt  was  to  increase  the  critical  trapping  velocity. 
However,  this  method  of  stirring  made  operation  quite  difficult.  It 
was  impossible  to  keep  the  distance  between  the  stirrer  and  the 
freezing  interface  close  and  constant.  If  the  distance  was  large, 
stirring  would  not  be  expected  to  influence  Vc>  If  the  stirrer  and 
the  interface  contacted  each  other,  either  the  stirrer  or  the  growth 
tube  was  damaged.  A  different  method  of  stirring  was  needed.  The 
method  of  horizontal  zone-refining  with  rotation  was  thus  adopted  to 
study  stirring  effect  on  rejection  of  particles. 

*  _ 


44-gauge  copper-constantan  thermocouple  with  0.15  mm  welded  junction 
diameter. 


The  method  of  horizontal  zone- ref ini ng  with  rotation  was 
developed  in  1966  by  Pfann,  Miller  and  Hunt  [42].  In  this  technique 
a  horizontal  tube  is  rotated  about  the  tube  axis.  If  particles  are 
present  in  the  molten  zone,  they  move  with  the  melt.  The  particles 
are  suspended  in  the  zone  and  are  not  always  in  contact  with  the 
freezing  interface.  A  Lepel  floating  zone  refiner  (model  FLZ-100) 
was  modified  to  suit  the  present  investigation.  The  molten  zone  was 
heated  with  a  loop  of  ni chrome  wire  (3  mm  width)  which  was  passed 
around  the  rotating  tube,  as  shown  in  Figure  8.  Two  heaters,  which 
could  be  controlled  independently,  were  needed  to  control  the  shape 
of  the  freezing  interface.  One  was  used  to  heat  the  zone;  the  other 
to  control  the  heat  flow  from  the  interface  into  the  growing  crystal. 
The  electrical  circuit  for  the  heating  system  was  similar  to  the  one 
used  in  the  Bridgman-Stockbarger  apparatus,  except  that  a  step-down 
transformer  (5.2  volts  c.p.  -  24A)  was  placed  between  the  constant 
voltage  transformer  and  the  auto- transformer. 


2.  Sample  Preparation 

Two  sizes  of  Pyrex  tube  (13  and  22  mm  O.D.)  with  30  cm  length 
were  used  to  study  the  effect  of  tube  diameter  on  particle  pushing 
with  stirring.  As  shown  in  Figure  9,  a  mixture  of  ul trasonically 
dispersed  particles  and  organic  melt  was  poured  into  a  tube  with  a 
slidable  Teflon  plug  about  3  cm  above  the  sealed  bottom.  Organic 
melt  without  particles  was  then  poured  on  top  of  this  mixture  which 


Melting  zone 
(15  mm) 


Figure  8  A  naphthalene  molten  zone  containing  carbon 
particles  traveling  in  the  right-hand 
direction  with  a  concave  freezing  interface 
(Only  left-hand  heater  used,  the  white 
material  appearing  in  the  melting  zone  is 
part  of  the  heater  support). 
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organic  compound 
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Mixture  of  particles 
and  organic 


j Teflon  plug 


Preparation  of  the  solid  charge  for. 
horizontal  zone-refining  with  rotation. 
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had  solidified.  The  purpose  of  the  slidable  plug  was  to  prevent  tube 
breakage  [42]  due  to  increases  in  volume  upon  melting. 

Because  the  initial  solid  charge  was  porous,  a  bubble  or  void 
appeared  at  the  top  of  the  molten  zone  during  horizontal  zone¬ 
refining  with  rotation.  As  the  zone  advanced,  the  porous  solid  was 
melted  and  a  nonporous  solid  formed  behind  the  zone.  The  void  thus 
gradually  increased  in  size  and  eventually  occupied  more  than  half  of 
the  zone  at  the  end  of  the  run.  The  void  adversely  affected  the 
particle  separation  process  when  it  occupied  more  than  one- third  of  the 
zone  volume.  Preliminary  in-situ  zone-melting  of  the  particle-free 
portion  of  the  solid  charge  (see  Figure  9)  was  therefore  used  prior  to 
the  particle-pushing  experiments  in  order  to  reduce  and  control  the 
size  of  the  void. 

3.  Experimental  Procedurp 

ror  a  typical  experiment,  a  charge  tube  prepared  as  above  was 
placed  horizontally,  rotated  about  its  own  axis,  and  moved  through 
a  stationary  heater.  The  particles  were  suspended  in  the  molten  zone 
and  moved  with  it  at  travel  rates  below  \JQ.  The  crystal  growth  rate 
equalled  the  movement  rate  under  steady-state  conditions.  The  growth 
rates  would  have  some  fluctuation  due  to  drafts  in  the  laboratory. 

The  rates  were  increased  step  by  step  (about  3  mm/hour  every  20 
minutes)  with  particle  trapping  again  observed  both  by  microscopic 
observations  and  by  sectioning  of  the  solidified  material.  The  errors 
in  the  measured  Vc  were  less  than  5%. 
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C.  Materials 
1 .  Particles 

Properties  of  the  particles  employed  for  these  experiments  are 
summarized  in  Table  2.  Carbon  particles  (bone  black)  were  extensively 
used  as  foreign  particles  in  this  study  because  they  appear  to  be  one 
of  the  common  impurities  i  organic  compounds.  Scanning  electron 
micrographs  of  these  particles  are  shown  in  Figure  10.  Special 
preparations  were  needed  to  take  these  pictures.  The  carbon  sample 
was  placed  on  a  piece  of  brass  which  was  coated  with  a  300-500  A  gold 
film.  It  appears  that  the  carbon  particles  were  actually  aoglomerates 
of  very  small  particles.  The  agglomerates  were  porous  and  of  irregular 
shape.  The  individual  particles  were  not  spherical  and  their  surfaces 
were  not  smooth.  These  agglomerates  were  directly  added  into  organic 
melts  for  early  Bridgman-Stockbarger  runs.  Later,  the  carbon  was  dis¬ 
persed  into  the  melt  using  ultrasonic  vibrations,  which  broke  the 
agglomerates  into  individual  particles  and  removed  most  of  the  gases 
entrapped  in  the  carbon. 

Copper  and  red  iron  oxide  were  also  studied  extensively  in  the 
naphthalene  melt  for  experiments  of  horizontal  zone-refining  with 
rotation.  Other  types  of  particles  such  as  silver,  zinc  and  silicon 
were  also  used  for  some  preliminary  studies. 

The  size  distributions  of  ultrasonically  dispersed  carbon,  iron 
oxide  and  copper  particles  were  determined  by  a  Mill i pore  ttMC  particle 
measurement  computer  system  [45].  Figures  11  and  12  show  the  distribu¬ 
tions  based  on  a  count.  The  average  diameters  were  1.0,  1.7  and 
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(A)  Carbon  particles  showing  agglomeration,  300X. 


I* -  200  A m  - 4 


(B)  Carbon  particle  enlarged  from  (A)  at  arrow,  6000X. 


Figure  10  Micrographs  of  carbon  particles. 


Size  distribution  for  carbon 


%  GREATER 

Figure  12  Size  distribution  for  copper  particles. 


26.0  um  for  iron  oxide,  carbon  and  copper,  respectively. 


2.  Organic  Compounds 

The  organic  compounds  studied  as  matrix  materials  were 
naphthalene,  camphor,  salol  and  benzophenone.  Their  properties  are 
summarized  in  Table  3.  Preliminary  results  indicated  that  only 
napnthalene  could  readily  push  these  particles  and  thus  it  was  used 
fcr  most  of  the  experiments.  Initially,  an  attempt  was  made  to  zone 
refine  naphthalene  and  camphor,  but  ironically  the  very  small  native 
particles  present  in  the  reagents  were  difficult  to  remove  thereby. 
Therefore,  these  chemicals  were  used  directly  in  experiments  without 
further  purification.  Some  of  these  native  particles  collected  from 
naphthalene  in  solidification  run  BNA-11B  are  shown  in  Figure  13. 
Fortunately,  such  particles  differed  sufficiently  from  those  added 
that  no  difficulty  was  encountered  in  measuring  Vc  of  added  particles. 
Details  of  zone-refining  camphor  and  naphthalene  are  discussed  in  the 
next  section,  along  with  gas  chromatographic  analyses  of  these  com¬ 
pounds  and  their  refined  products.  No  impurities  were  detected  in 
Baker  reagent  grade  naphthalene  by  conventional  gas  chromatography. 

D.  Solidified  Products 

In  addition  to  direct  observation  during  the  solidification 
experiment,  the  solidified  products  in  their  glass  containers  were 
examined  microscopically  (up  to  500  X)  both  through  the  glass  tube 
and  after  sectioning  to  determine  critical  velocities  and  particle 
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Figure  13  Native  particles  collected  from 
naphthalene  (Baker  analyzed 
reagent),  50X. 
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trapping  characteristics.  Samples  were  cut  either  traversely  or 
longitudinally  by  means  of  a  rotating  blade  which  is  mounted  on  a 
Micro-Matic  Precision  Watering  Machine  (Micromech  Mfg.  Corp.).  Some¬ 
times  the  naphthalene  slices  were  dissolved  in  acetone  and  the  solu¬ 
tion  filtered  through  a  1  urn  Mi  1 1 i pore  filter.  The  filter  plus 

collected  particles  was  placed  between  two  glass  slides  and  examined 
under  the  microscope. 

Zone-Refining  of  Organic  Compounds 

A  Fisher  zone  refiner  was  used  to  purify  organic  compounds  for 
the  particle-pushing  experiments.  Camphor  and  naphthalene  were  zone- 
refined  using  a  2  cm  O.D.  glass  tube  with  two  heaters  travelling 
upward  at  low  speeds  (2.5  to  10  mm/hour). 

1 .  Camphor 

Air  oxidation  of  camphor  occurred  during  the  initial  melting 
and  tube  filling  and  during  the  first  few  zone  passes.  The  brown- 
colored  oxidation  products  were  readily  removed  as  zone- refining 
continued.  Filling  the  tube  with  argon  did  not  completely  prevent 
the  discoloration,  because  some  air  was  entrapped  in  the  camphor. 

When  the  melt  was  uncolored,  the  freezing  interface  was  planar 
and  smooth  if  the  growth  rate  was  less  than  10  mm/hour.  Above  this 
rate  the  interface  became  dendritic.  Cells  were  observed  on  the 
interface  when  the  melt  was  colored  and  the  growth  rate  exceeded 
2.5  mi/hour.  Foreign  particles  originally  present  in  the  reagent  were 
very  difficult  to  remove.  Most  of  these  particles  were  at  first 
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pushed  by  the  planar  interface  at  9.7  mm/hour,  bet  were  then  subse¬ 
quently  incorporated  into  the  crystal  after  less  than  10  *  of 
pushing.  The  particle  separation  was  not  greatly  improved  when  the 
average  zoning  rate  was  decreased  to  2.5  mm/hour,  which  indicates  that 
there  was  a  fluctuating  freezing  rate  due  to  drafts  in  the  room,  power 
fluctuations,  and  the  intermittent  motion  of  the  heater. 

Analyses  by  gas  chromatography  indicated  that  the  Baker-grade 
DL-camphor  contained  about  five  different  trace  impurities.  Three  of 
these  impurities  were  removed  by  zone-refining  at  2.5  mm/hour  (RC-7). 
«e  identified  one  of  the  remaining  impurities  as  camphene,  the  con¬ 
tent  of  which  was  lower  in  the  refined  camphor.  The  concentration  of 
one  other  impurity  was  higher  in  the  product  than  in  the  starting 
material.  The  chromatograms  of  Baker,  refined  and  colored  camphor 
are  shown  in  Figure  14.  The  colored  oxidation  products  did  not  show 
up  in  the  chromatograms.  However,  they  coagulated  as  brownish  floes 
in  the  CS2  solution  used  for  chromatography,  after  the  solution  had 
stood  for  several  weeks.  When  reheated  in  air,  the  zone-refined 
camphor  became  yellow  while  no  change  occurred  if  oxygen  was  excluded, 
proving  that  air  oxidation  is  the  origin  of  the  colored  material. 


2.  Naphthalene 

This  compound  Is  much  more  stable  than  camphor  when  heated  in 
At  a  zone  speed  of  4  mm/hour  (RN-6),  the  native  particles 
present  in  Baker-analyzed  reagent  were  difficult  to  remove.  The 
freezing  interface  was  smooth  and  planar.  Occasionally,  surface- 


(C)  Colored  portion  (105.4  mg/100  cc  CS2). 


Figure  14  Chromatograms  of  camphor  (Aerograph  HY-F1 , 
Model  600,  column  5  ft  x  l/8  in,  5#  SE-30 
60/80  A/W  DMCS  Chrome  W,  each  injection 
4.5  Al) . 
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driven  flows  caused  particles  to  vigorously  circulate  around  bubbles. 
Gas  chromatographic  analyses  revealed  no  impurities  in  the  naphthalene 
(Figure  15),  even  with  the  use  of  a  boosting  technique  and  high 
temperature  (360°C)  injection.  Thus,  the  impurity  content  must  be 
less  than  0.1%,  which  is  the  minimum  detectable  value  for  this  gas 
chromatograph.  The  zone-refined  (bottom)  and  slightly  yellow  (top) 
portions  (RN-9)  gave  chromatograms  similar  to  the  reagent,  with  no 
impurity  peak.  When  0.5  wt  %  of  anthracene,  a  likely  impurity,  was 
added  to  the  naphtialene,  the  peak  showed  up  as  shown  in  Figure  15(B). 


Figure  15  Chromatograms  of  naphthalene 
(Aerograph  HY-F1 ,  Model  600D, 
same  column  as  used  in  camphor). 


CHAPTER  IV 

RESULTS  AND  DISCUSSIONS 


A  systematic  experimental  study  was  originally  planned  to 
determine  the  effect  of  two  important  parameters-interfacial  tempera¬ 
ture  gradient  and  stirring-on  the  critical  trapping  velocity.  Many 
preliminary  experiments  were  conducted  to  find  a  suitable  particle- 
organic  matrix  pair.  The  ideal  pair  should  have  a  measurable  critical 
velocity,  non-faceted  solid-liquid  interface  and  easy  handling  for 

sample  preparations.  The  pair  of  carbon  and  naphthalene  was  found  to 
fulfill  most  of  the  requirements. 

Many  phenomena,  which  help  in  the  understanding  of  particle 
pushing,  were  discovered  during  the  experiments  and  later  treatment  of 
the  solidified  products.  Experiments  were  designed  not  only  to  pro¬ 
vide  quantitative  data  on  the  effect  of  temperature  gradient  and 
stirring,  but  also  to  assist  in  the  discovery  of  new  and  important 
phenomena  through  direct  microscopic  observation. 

The  experimental  results  obtained  from  the  vertical  Bridgman 
growth  method  are  discussed  first.  Those  obtained  from  horizontal 
zone- refining  with  rotation  follow.  A  summary  of  results  is  given  at 
the  end  of  the  discussion  for  each  experimental  method. 


A-  Bn.dgman-Stockbarqer  Technique 


1  •  General 

Consider  inert  particles  more  dense  than  the  organic  melt  added 
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into  the  upper  portion  of  the  growth  tube  (see  Figure  6).  The 
particles  are  acted  upon  by  various  forces  (such  as  gravity,  viscous 
drag  and  bouyancy)  and  settle  in  the  melt.  Most  of  them  arrive  on 
the  horizontal  solid-liquid  interface,  and  a  new  steady-state  inter¬ 
face  in  contact  with  the  particles  is  formed.  Very  fine  particles 
are  suspended  in  the  melt  due  to  Brownian  motion  and  gentle  convective 
currents  which  exist  in  the  melt  about  1  cm  above  the  freezir.g  inter¬ 
face.  As  solidification  starts  at  slow  speeds,  the  particles  covering 
the  interface  are  swept  along  as  the  freezing  interface  advances. 

The  moving  interface  also  contacts  the  fine  particles  suspended  in 
the  melt  and  pushes  them.  As  the  growth  rate  is  gradually  increased, 
a  speed  is  reached  at  which  some  particles  cease  to  be  pushed  by  the 
interface  and  are  trapped  in  the  growing  crystal.  With  many  different 
sizes  and  shapes  of  particles,  the  observed  transition  from  pushing  to 
trapping  takes  place  over  a  wide  range  of  velocities. 

The  significant  experiments  performed  in  the  vertical  Bridgman 
method  are  tabulated  in  Table  4  (salol)  and  Table  5  (naphthalene), 
together  with  the  observed  results. 

2.  Solid-Liquid  Interface  Morphology 

The  solid-liquid  interface  was  maintained  by  heating  the  uDper 
portion  of  the  growth  tube  and  cooling  the  lower  portion  of  the  tube. 
The  shape  of  this  interface  was  significantly  affected  by  power  input 
to  the  heater  and  by  coolant  temperature.  Three  kinds  of  interface 
shapes  were  observed  in  the  vertical  Bridgman  growth  method,  as  shown 


TABLE  4 

SUMMARY  OF  VERTICAL  BRIDGMAN  EXPERIMENTS  FOR  CARBON-SALOL  SYSTEM^1 ^ 
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in  Figure  16.  At  a  given  growth  rate,  the  freezing  interface  moved 
upward  into  the  heating  zone  when  the  coolant  temperature  was  lowered 
and  became  convex  (Figure  16A).  On  the  other  hand,  when  the  heater 
power  was  increased,  the  freezing  interface  moved  downward  into  the 
cooling  zone  and  became  concave  (Figure  16C).  A  planar  interface 
could  be  obtained  by  adjusting  conditions  so  that  the  interface  was 
in  the  vicinity  of  the  junction  of  the  heater  and  cooler  (Figure , 16B) . 

The  shape  of  the  interface  influences  the  critical  trapping 
velocity.  A  convex  interface  caused  most  of  the  particles  to  fall 
into  the  groove  created  between  the  interface  and  glass  wall.  The 
particles  at  the  glass-solid-liquid  interface  were  trapped  sooner 
than  those  in  the  center  of  the  tube.  The  former  gave  a  Vc  of  18  mm/ 
hour,  and  the  latter  a  value  of  21  mm/hour  for  ul trasonical ly  dis¬ 
persed  carbon  in  naphthalene.  Figure  17  shows  carbon  trapped  between 
the  crystal  and  the  tube  wall. 

The  freezing  interface  of  naphthalene  was  always  smooth  under 
the  imposed  temperature  gradients  of  this  study,  even  though  naphtha¬ 
lene  has  a  large  entropy  of  fusion  (ASf  =  53.4  J/mole  K).  The 
interface  morphology  of  salol  was  influenced  by  the  growth  rate  in 
addition  to  the  temperature  gradient.  At  high  temperature  gradients 
the  interface  became  smooth  when  the  growth  rate  was  less  than  6  mm/ 
hour,  but  it  was  faceted  when  the  growth  rate  increased  above  8  mm/ 
hour. 

3.  Temperature  Gradients 

Direct  addition  of  carbon  into  the  melt  of  naphthalene  was  used 
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to  study  the  effect  of  temperature  gradient  on  the  critical  velocity. 
With  varied  gradients,  the  critical  velocities  were  all  about 
12  mm/hour  as  shown  in  Table  6.  Higher  gradients  were  obtained  by 
means  of  decreasing  the  temperature  of  the  cooling  fluid  and 
increasing  the  heater  power  input.  The  measured  interfacial  gradients 
are  only  approximate. 

At  the  end  of  this  preliminary  study,  Cisse  and  Bolling  [9] 
presented  similar  results  based  on  particle-water  systems.  Their 
results  also  indicated  that  the  critical  velocity  does  not  depend 
significantly  on  the  interfacial  temperature  gradient.  Because  of 
the  foregoing  evidence  and  similar  conclusions  based  on  the  organic 
system,  we  did  not  pursue  this  subject  further. 

Theoretical  studies  of  the  factors  influencing  the  interfacial 
temperature  gradient  were  also  made  during  this  research.  A  paper 
was  published  based  on  part  of  this  study  and  is  attached  as 
Appendix  A.  The  influence  of  the  length  and  radius  of  a  crystal  on 
the  interfacial  temperature  gradient  was  estimated  by  use  of  a  one¬ 
dimensional  heat  transfer  calculation.  It  was  found  that  the  teinpera- 
t'.re  gradient  continues  to  change  for  quite  long  crystals  when  the 
Biot  number  H  =  hR/k  is  small  (where  h  is  the  heat  transfer  coef¬ 
ficient  from  the  crystal  surface,  R  is  crystal  radius,  and  k  is  the 
thermal  conductivity),  especially  for  good  conductors,  but  does  not 
depend  on  crystal  length  significantly  when  the  Biot  number  is  large. 
For  our  Bridgman  growth  of  naphthalene,  H  was  about  1  and  the  critical 
length  was  then  equal  to  R.  Therefore,  a  length  of  tube  equal  to  its 
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TABLE  6 

EFFECT  OF  TEMPERATURE  GRADIENT  ON  CRITICAL  VELOCITY 
{Naphthalene-Carbon  System)^  ^ 


Run  No. 

Bath 

Temperature, 

°C 

Heater 

Input, 

watts 

Temperature 

Gradient, 

°C/cm(2) 

Cri tical 
Veloci ty , 
mm/hr. 

BN-5 

-20 

6.8 

~  27 

~1 2 

BN- 3 

0 

5.0 

~  6 

~1 2 

BN- 7 

10 

4.5 

- 

~1 2 

BN-8 

20 

4.5 

- 

~1 2 

^Carbon  added  directly  into  the  melt. 

^Measured  interfacial  temperature  gradient  in  the  melt  at  the 
interface. 
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diameter  was  submerged  in  the  cooler  at  the  beginning  of  each  experi¬ 
ment,  so  that  constant  temperature  gradients  were  maintained  at  the 
freezing  interface  during  controlled  solidification. 


4.  S 1 1 rri nq 

It  was  hoped  that  stirring  of  the  melt  would  increase  the 
critical  velocity.  Stirring  should  enhance  the  supply  of  melt  to  the 
contact  region  between  the  particles  and  the  interface,  move  particles 
along  the  interface  giving  them  less  time  to  be  trapped,  and  keep 
many  more  of  the  particles  in  suspension  and  away  from  the  interface. 
Because  of  its  very  low  vapor  pressure,  salol  was  originally  used  for 
the  stirring  experiments  with  carbon.  The  experiment  without  stirring 
gave  Vc  -  7  mn,/hour.  however,  the  runs  with  agitation  did  not  give 
satisfactory  results.  When  carbon  was  added  directly  into  the  melt 
(BSA-9/1),  most  of  the  carbon  was  mediately  incorporated  into  the 
growing  salol  and  made  a  dark  band  as  shown  in  Figure  10  (three 
additions  were  made  in  this  run),  when  the  carbon  was  dispersed  into 
the  melt  ultrasonically  (USA-12A),  the  interface  trapped  particles 
quite  uniformly,  resulting  in  a  nearly  uniformly  dark  solid.  Other 
kinds  of  particles,  Fe.O,  (BSA-13-1),  Si  (BSA-,4)  and  Zn  (BSA-13-2), 
were  also  used  which  were  reported  to  be  easily  separated  from  salol 
By  Uhlmann  et  al.  [6],  However,  we  did  not  obtain  the  expected 
results  even  without  stirring  (see  Table  4). 

Therefore,  the  naphthalene-carbon  system  was  used  henceforth 
in  the  study  of  stirring,  due  to  the  ready  separation.  Inserting  a 
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stirrer  in  the  growth  tube  made  operation  quite  difficult.  It  was 

impossible  to  keep  the  distance  between  the  stirrer  and  the  interface 

constant.  If  the  distance  was  large,  stirring  would  not  be  expected 

to  influence  V  .  If  the  stirrer  and  the  interface  contacted  each 
c 

other,  either  the  growth  tube  or  the  stirring  rod  were  damaged. 
Although  we  used  a  stepless  jack  to  adjust  the  position  of  the 
stirrer  frequently  during  changes  in  the  growth  rate,  this  distance 
was  very  difficult  to  control  and  made  the  effect  of  stirring 
uncertain.  However,  one  stirring  run  (BNA-8A)  gave  Vc  =  30  mm/hour, 
which  was  quite  encouraging. 

Another  phenomenon  was  observed  as  follows.  When  stirring  was 
initiated,  the  shape  and  position  of  a  steady  salol  interface  changed 
drastically.  The  planar  interface  became  convex  into  the  solid  due 
to  higher  heat  transfer  rates  through  the  glass  wall.  Particles  moved 
up  the  bulge  and  circulated  into  the  melt  (motion  pictures  were  taken 
of  this  aspect).  After  the  interface  reached  steady-state,  it  was 
concave  into  the  solid  making  direct  microscopic  observation  of 
particle  pushing  difficult.  The  melt  was  decanted  (BSA-6A,  7A)  and 
the  interface  was  observed  thereby  to  consist  of  many  large  facets, 
as  shown  in  Figure  19.  Some  grains  protruded  up  to  1.5  mm.  The 
deepest  point  in  the  concave  interface  was  about  9  mm  from  the 
periphery. 

5.  Matrix-Particle  Dependence 

Previous  data  [6]  indicated  that  there  were  variations  in  V 

c 
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among  the  particle  materials  in  a  given  matrix  material,  and  among 
the  matrix  materials  for  a  given  particle.  We  used  a  different, 
improved  experimental  method  and  observed  that  carbon  particles  were 
readily  separated  from  naphthalene,  but  were  very  difficult  to  remove 
from  camphor  (Vc  less  than  2  mm/hour).  A  summary  of  results  is  listed 
in  Table  7. 

6 .  Gas  Bubbles  and  Surf ace-Driven  Flows 

The  solubility  of  gases  in  solids  is  normally  much  lower  than 
in  the  corresponding  melts.  Segregation  takes  place  at  the  freezing 
interface  and  increases  the  gas  concentration  next  to  the  interface. 

As  a  result,  gas  bubbles  often  nucleate  and  grow  on  the  interface 
during  solidification  [49].  Because  temperature  varies  in  the  melt, 
the  surface  tension  varies  along  the  gas-liquid  interface  and  can 
produce  vigorous  circulation.  In  the  present  experiments,  bubbles 
gave  rise  to  vigorous  flows  which  circulated  the  particles  around  the 
bubbles.  Figure  20  shows  vigorous  circulation  of  carbon  particles 
around  a  bubble  growing  on  the  interface  of  naphthalene  (frame  from 
motion  picture). 

These  gas  bubbles  had  a  large  influence  on  the  separation  of 
foreign  particles  from  the  matrix  materials.  Gas  "worms"  were  often 
created  in  the  solid  behind  the  bubbles  on  the  interface,  as  shown  in 
Figure  21.  These  bubbles  sometimes  floated  away.  Sometimes  a  bubble 
moved  rapidly  back  into  its  gas  worm  or  into  a  sudden  crack  in  the 
solid  and  deposited  a  large  number  of  particles  in  the  solid.  In  such 
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TABLE  7 


CRITICAL 

TRAPPING  VELOCITIES  IN  VERTICAL 
BRIDGMAN  CRYSTAL  GROWTH 

Organic  Matrix 

Particle 

vc 

mm/ hour 

Remarks 

Naphthalene 

C 

21 

U1  trasonical ly 
dispersed. 

16 

Trapping  next  to 
glass  wall. 

12 

Agglomerate 

particles. 

Cu 

18 

14 

Trapping  next  to 
glass  wall. 

Salol 

C 

7 

^e2^3 

not  pushed 
at  6 

Zn 

not  pushed 
at  6 

Si 

not  pushed 
at  6 

Camphor 

C 

not  pushed 
at  2 
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cases  it  was  very  difficult  to  know  whether  the  particle  incorpora¬ 
tion  was  due  to  a  true  critical  velocity  or  to  the  influence  of 
bubbles. 

Some  interesting  observations,  not  related  to  particle  pushing, 
were  also  made  on  bubble  nucleation  and  growth.  In  the  carbon- 
naphthalene  system,  bubbles  occurred  at  a  change  of  crystal  growth 
rate,  or  when  the  growth  rate  exceeded  16  mm/hour  (BNA-4).  When  the 
growth  tube  was  evacuated  at  0.3  Torr,  bubbles  did  not  appear  unitl 
V  =  26  mm/hour  (BNA-6V,  7 V).  However,  no  bubbles  formed  when 
copper  foreign  particles  were  employed  in  naphthalene  with  a  0.3 
Torr  evacuation  even  at  V  =  32  mm/hour  (BNA-9V  Cu).  Hence,  nuclea 
tion  of  bubbles  depended  on  the  type  of  particles  present  and  the 

amount  of  gas  in  the  matrix  material. 

Stirring  the  melt  reduced  bubble  formation  (BNA-1A,  2A,  3A), 
but  did  not  eliminate  the  problem  (BNA-8A) .  Bubbles  also  formed  on 
the  freezing  interface  of  salol  under  some  conditions. 

A  theoretical  treatment  of  bubble  nucleation  during  solidifica¬ 
tion  is  presented  in  Chapter  VI. 

7.  Grain  Surfaces 

The  solidified  material  in  this  study  normally  consisted  of  2 
to  4  grains  for  naphthalene  and  about  12  grains  for  salol  when  using 
a  13  mm  O.D.  glass  tube.  Several  single  grains  were  separated  from 
one  another  and  examined  under  a  microscope.  In  salol,  grain 
boundaries  contained  many  steps,  cavities,  gas  worms  and  preferentially 
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trapppd  particles,  as  shown  in  Figure  22.  However,  the  surfaces  of 
naphthalene  grains  primarily  exhibited  steps  without  preferential 
trapping  of  particles  (Figure  23).  The  latter  appears  to  confirm  the 
results  of  Cisse  and  Bolling  [9],  in  which  they  found  particles  on 
grain  boundary  grooves  are  pushed  to  higher  velocities  than  those  at 
flat  grain  surfaces  of  an  ice  water  interface. 

8.  Growth  Orientation 

Examination  of  the  solidified  products  of  many  runs  indicated 
that  one  side  of  the  interface  began  trapping  particles  much  sooner 
than  the  other  side,  as  shown  in  Figure  24A.  This  leads  one  to 
believe  that  certain  grain  orientations  may  push  particles,  while 
other  grain  orientations  incorporate  particles  at  the  same  growth  rate. 
The  transverse  cut  of  the  solidified  ingot  (BNA-8A)  showed  that  one 
grain  trapped  particles,  while  the  other  grains  rejected  particles 
(see  Figure  24B). 

Incidentally,  the  rejection  of  soluble  impurities  was  also 
observed  to  depend  on  growth  orientation.  In  zone-refined  camphor 
one  grain  was  light  yellow,  while  the  adjacent  grain  was  colorless. 

9.  Surfaces  of  Solidified  Products 

Like  the  surfaces  of  the  grain  boundaries,  the  cylindrical 
surfaces  where  the  solid  contacted  the  glass  tube  were  not  always 
smooth.  Frequently,  they  were  rough  and  contained  steps,  bumps,  tiny 
crystals  (Figure  25),  and  gas  bubbles  and  worms  (Figure  26).  As 
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Figure  25  Surfaces  of  solidified  naphthalene 
at  tube  wall. 
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mentioned  before,  carbon  particles  were  trapped  between  these  surfaces 
and  the  glass  wall  at  a  lower  critical  velocity  than  at  the  advancing 
grain  surfaces,  especially  with  convex  interfaces.  This  indicates 
that  the  results  of  Uhlmann,  et  al.  [6]  obtained  with  thin  layers  of 
melt  held  between  microscope  slides  might  have  been  strongly  influ¬ 
enced  by  the  glass  slides  and  might  not  have  been  characteristic  of  a 

free  interface. 

10.  Cracking  and  Healing 

Because  the  coefficient  of  thermal  expansion  differs  for  a 
solid  and  its  growth  container,  thermal  stresses  form  during  the 
period  of  cooling  after  solidification.  Sometimes  thermal  stresses 
were  sufficient  to  cause  cracking  of  solidified  material  both  in  the 
horizontal  and  vertical  directions.  The  cracks  did  not  reach  the 
freezing  interface.  Therefore,  particles  were  not  preferentially 
trapped  at  these  locations.  A  vertical  crack  was  sometimes  difficult 
to  distinguish  from  a  grain  boundary.  Occasionally  a  horizontal 
crack  healed  to  form  a  veil  of  microscopic  gas  bubbles,  as  shown  in 
Figure  27.  Similar  observations  were  made  in  movement  of  inclusions 

in  KI  crystals  [50]. 

11 .  Bouncing  Particles 

An  interesting  phenomenon  was  observed  during  zone-refining  of 
camphor.  Some  of  the  small  foreign  particles  originally  present  in 
the  reagent  bounced  back  into  the  melt  as  soon  as  they  touched  the 


Veil  of  microscopic 
gas  bubbles  in  healed 
portion 


Unhealed 

portion 


Crack  (black  area)  and  partially  healed 
crack  (a  veil  of  microscopic  gas  bubbles) 
in  solidified  naphthalene  (50X,  NBA-4). 
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freezing  interface.  Others  bounced  into  the  melt  following  a  rapid 
approach  to  the  interface  without  actually  touching  it.  This  is 
suggestive  of  electrostatic  effects  or  of  irregular  convection 
currents.  Further  study  is  warranted. 

12.  Summary 

In  addition  to  previously  established  parameters  such  as 

particle  size  and  shape,  particle  roughness,  and  grain  boundaries, 

we  have  found  that  bubbles,  interface  shape,  growth  container  and 

the  growth  orientation  all  significantly  influence  the  critical 

velocity  in  the  carbon-naphthalene  system.  Because  of  this  complex 

dependence,  it  was  never  observed  that  all  particles  were  clearly 

incorporated  at  one  location  in  the  solidified  products.  In  fact, 

the  particles  were  trapped  over  a  finite  distance.  The  measured 

values  of  V  varied  from  21  to  35  mm/hour  (BNA-7V)  without  stirring 
c 

the  melt.  The  growth  of  bubbles  at  the  freezing  interface  was 
eventually  prevented  by  evacuating  the  growth  container.  The  freezing 
interface  had  to  be  kept  slightly  concave  into  the  growing  solid  to 
prevent  particles  from  being  trapped  between  the  crystal  surface  and 
the  container.  The  interfacial  temperature  gradient  did  not  affect 
the  critical  trapping  velocity  wi thin  the  range  investigated.  The 
Vc  for  trapping  of  carbon  by  solidifying  salol  was  7  mm/hour,  but 
particles  of  Zn,  Si  and  Fe203  were  not  readily  pushed.  These  results 
are  contrary  to  data  by  others  [6]  in  which  a  horizontal  thin  film 
held  between  two  glass  slides  with  a  vertical  freezing  interface  was 
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employed.  Stirring  might  have  increased  Vc,  but  a  new  method  to 
perform  stirring  was  needed  to  replace  the  unsuccessful  technioue  of 
inserting  a  stirrer. 

B.  Horizontal  Zone-Refining  with  Rotation 
1 .  Interface  Shape 

Tube  rotation  averaged  out  radial  asymmetries.  The  shapes  of 
the  freezing  and  melting  interfaces  in  a  horizontal  rotating  zone 
depended  on  the  growth  rate.  The  variations  of  interface  shapes  are 
shown  in  Figure  ?8.  Planar  interfaces  were  obtained  only  at  zero 
or  very  slow  growth  rates.  When  the  growth  rate  was  increased,  the 
latent  heat  liberated  at  the  freezing  int^-face  had  to  be  conducted 
through  the  crystal.  Since  cooling  took  place  only  at  the  rotating 
tube  wall,  the  freezing  interface  gradually  became  more  and  more  con¬ 
cave.  Meanwhile,  the  melting  interface  moved  closer  to  the  heater 
(heat  power  kept  constant)  and  became  increasingly  convex,  likewise 
because  of  latent  heat  consumption.  At  high  growth  rates  (  >90  mm/ 
hour),  insufficient  time  was  available  to  completely  melt  through  the 
feed  and  the  melting  interface  began  contacting  the  freezing  interface, 
as  illustrated  in  Figure  28D.  Because  of  this,  the  critical  velocities 
of  certain  particle-organic  systems  with  high  Vc  were  difficult  to 
determine  even  when  heating  power  was  continuously  increased  during 
the  experiment  (for  example,  carbon-naphthalene  in  a  19  mm  I.U.  growth 
tube,  HRL  -  5  Cu). 

As  the  freezing  rate  of  naphthalene  was  increased,  small 
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Figure  28  Change  of  zone  shape  fith  trave J  jjjj  Jn 
horizontal  zone-refining  of  ^Phthalene 
with  rotation  (50  RPM.  19  mm  I.D.  tube). 
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amounts  of  carbon  particles  gradually  accumulated  on  the  central 
portion  of  the  concave  interface  and  adhered  there.  These  particles 
were  then  trapped  into  the  growing  crystal  when  the  growth  rate 
reached  \IQ.  The  length  of  crystal  over  which  particles  were  entrapped 
depended  on  the  growth  rate  after  exceeding  V  .  The  trapping  length 
was  shorter  as  the  growth  rate  was  increased.  Periodic  trapping 
was  observed  over  a  finite  length  when  the  growth  rate  was  kept  at  V  . 
In  other  words,  accumulation  of  particles  on  the  concave  interface 
and  trapping  of  particles  into  the  crystal  took  place  alternatively 
at  until  all  particles  in  the  rotating  zone  were  incorporated. 

Figure  29  shows  that  carbon  particles  were  pushed  from  the  left  end 

to  the  right  end  of  a  naphthalene  ingot  by  the  rotating  zone-refining 
method. 

Red  iron  oxide  particles  gave  trapping  characteristics  similar 
to  carbon.  Both  were  entrapped  along  the  axis  of  growing  naphthalene. 
However,  the  pushing  and  trapping  phenomena  were  different  for  copper 
particles.  At  Vc  the  particles  did  not  accumulate  at  the  central 
part  of  the  interface,  but  were  incorporated  into  the  crystal  rather 
uniformly. 

The  effect  of  interface  shape  on  particle  pushing  was  also 
studied.  When  the  freezing  interface  was  kept  planar  by  means  of  a 
second  heater,  the  results  were  slightly  different  from  those  obtained 
for  a  concave  interface.  Carbon  particles  were  trapped  uniformly  into 
the  growing  naphthalene  and  Vc  was  increased  from  30  to  35  mm/hour 
(HRT-6,  7  and  8).  A  planar  interface  was  not  easily  obtained  because 
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igure  29  Horizontal  zone-refined  naphthalene 
with  carbon  particle q  pushed  from  the 

aleJV°  *h?  right  end  (vc=30  mm/hr 
and  5  RPM  in  a.  13  mm  OD  tube). 
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each  heater  required  a  different  heat  input  dependent  on  growth  rate 
and  rotation. 

2.  Rotation 

A  minimum  speed  of  rotation  was  first  established  in  order  to 
keep  most  of  the  particles  suspended  in  the  rotating  zone  for  each 
particle-organic  matrix  pair.  The  rotation  was  then  increased  to 
determine  the  effect  of  rotation  speed  on  \IQ.  Tables  8,  9  and  10  are 
summaries  of  results  for  carbon,  copper  and  red  iron  oxide  in  a 
naphthalene  matrix. 

Two  kinds  of  Vc  are  shown  in  these  tables.  One  is  the  V  for 
particle  trapping  into  the  solid,  the  other  is  for  trapping  between 
the  tube  wall  and  the  solid.  Copper  particles  gave  a  strong 
dependence  of  Vc  on  rotation  speed.  Increase  of  rotation  increased 
the  critical  trapping  rate  of  copper  into  the  solid.  On  the  other 
hand,  increase  of  rotation  increased  the  centrifugal  force  causing 
some  of  the  larger  particles  to  stick  to  the  tube  wall.  These 
particles  were  then  engulfed  next  to  the  tube  wall.  Figures  30,  31 
and  32  indicate  the  variations  of  Vc  in  the  two  positions  with  rotation 
speed  for  copper,  iron  oxide  and  carbon.  For  example,  in  the  copper- 
naphthalene  system  Vc  in  the  crystal  increased  with  rotation,  while 
Vc  next  to  the  wall  decreased  with  rotation.  The  intersection  of 
these  two  curves  gives  an  optimum  rotation  speed  under  which  a 
maximum  V£  at  the  crystal  is  reached  with  no  trapping  next  to  the  wall. 

The  optimum  rotation  speeds  and  corresponding  Vc  are  tabulated  in 
Table  11. 
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TABLE  8 

EFFECT  OF  ROTATION  RATE  ON  CRITICAL  FREEZING  RATE  V 
FOR  TRAPPING  OF  CARBON  IN  NAPHTHALENE 
USING  HORIZONTAL  ZONE  REFINING 


Run 

Tube  Size 
mm  I.D. 

Rotation 

RPM 

Cri ti cal 

Crystal 

Velocity,  mm/hour 

Wall 

HR- 6 

10.5 

1.5 

29 

30 

HR-8 

10.5 

1.4 

28 

30 

HR- 2 

10.5 

A 

34 

no 

HR- 5 

10.5 

4 

30 

no 

HR-7 

10.5 

16 

29 

30 

HR-3 

10.5 

33 

32 

32 

HR-4 

10.5 

33 

30 

31 

HR-15M 

10.5 

45 

33 

28 

HR-16M 

10.5 

45 

34 

29 

HR- 23 

10.5 

80 

30 

20 

URL- 1 

19 

10 

24 

28 

HRL-2 

19 

30 

36 

36 

HRL-3 

19 

60 

37 

trapped 
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TABLE  9 

EFFECT  OF  ROTATION  RATE  ON  THE  CRITICAL  FREEZING  RATE  V 
FOR  TRAPPING  OF  COPPER  IN  NAPHTHALENE  c 

USING  HORIZONTAL  ZONE  REFINING 


Run 

Tube  Size 
mm  I.D. 

Rotation 

RPM 

Critical  Velocity,  mm/hour 

Crystal 

Wall 

HR- 11 

10.5 

44 

47 

no 

HR- 12 

10.5 

44 

51 

no 

HR-13 

10.5 

44 

48 

no 

HR- 19 

10.5 

50 

66 

80 

HR-18 

10.5 

58 

70 

68 

HR- 17 

10.5 

70 

77 

65 

HR- 14 

10.5 

100 

82 

61 

HRL-5 

19 

44 

<  90 

53 

HRL-8 

19 

50 

<  90 

90 
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Run 

HR-21  X 
HR-20X 
HR-22X 
HRL-7X 
HRL-6X 


TABLE  10 

EFFECT  OF  ROTATION  ON  Vc  FOR  IRON  OXIDE-NAPHTHALENE 
IN  HORIZONTAL  ZONE  REFINING 


Tube  Size 

Rotation 

Critical 

Velocity,  mm/hour 

mm  I.D. 

RPM 

Crystal 

Wall 

10.5 

8 

28 

33 

10.5 

40 

35 

34 

10.5 

80 

36 

trapped 

19 

8 

32 

no 

19 

40 

46 

20 

40  50  60  70  80  90  100 

Rotation  Speed,  RPM 


Figure  30  The  influence  of  tube  rotation  on 
incorporation  of  copper  particles 
by  naphthalene  during  horizontal 
zone  melting  (10  mm  I.D.  tube). 


Rotation  Speed,  RPM 


Figure  31  The  influence  of  tube  rotation  on 

incorporation  of  iron  oxide  particles 
by  naphthalene  during  horizontal  zone 
melting  (10  mm  I.D.  tube). 


TABLE  11 


OPTIMUM  ROTATION  RATES  FOR  SEPARATION  OF  PARTICLES 
FROM  NAPHTHALENE  BY  A  HORIZONTAL  ZONE  REFINER 
WITH  ROTATION 

Tube  Size  Optimum  Rotation 

Particles  mm  I.D.  RPM _ mm/hour 


3.  Tube  Diameter 

Tables  8,  9,  10  and  11  also  give  the  experimental  results 
obtained  for  a  larger  growth  tube  (19  mm  I.D.,  22  mm  O.D.).  By 
increasing  the  tube  diameter  by  a  factor  of  two,  Vc  was  increased 
about  10%  for  iron  oxide,  20%  for  carbon  and  more  than  30%  for 
copper  under  the  respective  optimum  rotation  rates.  Furthermore,  the 
optimum  rotation  rates  for  the  large  tube  were  reduced  at  least  10% 
from  that  of  the  small  tube. 

4.  Bubbles  and  Void 

The  dissolved  gases  in  the  organic  compound  are  continuously 
rejected  into  the  melt  during  solidification.  These  gases  grew  into 
bubbles  at  the  advancing  interface  in  the  vertical  Bridgman 
experiments.  In  horizontal  zone  refining  with  rotation,  however,  the 
gases  were  transported  into  the  melt  sufficiently  to  avoid  nucleation 
and  growth  at  the  interface.  These  gases  were  liberated  in  a  gas 
space  or  void  located  at  the  top  of  the  rotating  zone.  The  void 
helped  in  observation  of  particles  at  the  interface.  Carbon  and  iron 
oxide  particles  always  adhered  at  the  center  of  the  interface  before 
they  were  incorporated  into  the  crystal.  Sometimes,  adhesion  of 
carbon  particles  on  the  interface  gave  a  spiral  pattern. 

Bubbles  appeared  at  the  interface  only  at  slow  rotation  rates 
and  high  growth  rate  during  separation  of  carbon  from  naphthalene. 
They  did  not  give  difficulties  in  determination  of  Vc  because  they 
be^an  occurring  at  about  30  mm/hour.  Bubbles  were  not  encountered  at 


the  interface  in  any  other  horizontal  zone  melting  experiments. 


5.  Summary 

Horizontal  zone  refining  with  rotation  achieved  the  goal  of 
increasing  the  critical  trapping  velocity  of  particles.  Sufficient 
data  were  gathered  to  permit  design  of  a  separation  process  and  to 
make  laboratory-scale  production  of  particle-free  organics.  This  new 
separation  process  is  simple  and  effective  and  gives  significant 
improvements  and  advantages  over  the  vertical  Bridgman  method.  The 
advantages  are  listed  as  follows:  (1)  Vc  is  increased  about  50%  for 
carbon  particles  and  at  least  300%  for  copper  in  naphthalene;  (2)  clear- 
cut  separation  is  possible;  (3)  trapping  between  the  crystal  and  the 
tube  wall  can  be  prevented;  (4)  bubbles  do  not  grow  on  the  freezing 
interface,  so  that  the  growth  tube  does  not  have  to  be  evacuated;  and 
(5)  operation  is  simple,  effective  and  economical.  The  economics  of 
this  separation  process  will  be  discussed  in  the  next  chapter. 

An  optimum  rotation  speed  was  determined  for  each  particle- 
naphthalene  system,  above  which  a  considerable  amount  of  particles  were 
trapped  next  to  the  glass  wall  and  below  which  trapping  was  preferred 
in  the  crystal.  Further  discussion  on  horizontal  rotation  is  given  in 
Chapter  VI. 


CHAPTER  V 

f^ENTIAL  APPLICATIONS  AND  ECONOMICS 


In  addition  to  investigation  of  particle  pushing  and  trapping 
by  a  freezing  interface,  the  final  goals  of  this  work  were  to  make 
particle- free  organic  chamicals,  especially  those  which  are  solid  at 
room  temperature.  In  addition,  several  potential  applications  were 
discovered  during  the  course  of  this  work.  Some  of  them  were  further 
studied  experimentally.  The  results  and  economics  of  the  developed 
process  are  given  as  follows. 

A.  Separation  of  Mixed  Particles 

The  results  shown  in  Chapter  IV  and  previous  data  [6]  have 

demonstrated  that  each  kind  of  particle  has  a  different  freezing  rate 

Vc  in  a  given  matrix  material.  In  other  words,  a  mixture  of  different 

kinds  of  particles  may  be  separated  in  a  properly  selected  organic 

compound  if  an  increasing  growth  rate  is  employe'  during  controlled 

solidification.  The  results  of  horizontal  zone  refining  with 

rotation  gave  V£  =  30  mm/hour  for  carbon  and  Vc  =  70  mm/hour  at  55  RPM 

for  copper  in  naphthalene.  Because  of  their  wide  difference  in  V  ,  it 

c 

should  be  possible  to  separate  a  mixture  of  carbon  and  copper  from 
one  another  by  freezing  of  naphthalene.  Run  HR-15M  accomplished  this 
goal.  When  the  freezing  rate  was  progressively  increased  from  25 
to  75  mm/hour  at  44  RPM  in  a  period  of  4  hours,  carbon  particles 
started  to  be  trapped  at  34  mm/hour.  All  of  the  carbon  was  incorporated 
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into  the  central  portion  of  the  crystal  over  a  length  of  2  cm. 

Copper  particles  were  then  trapped  uniformly  at  54  mm/hour,  until 
the  end  of  the  run  (75  mm/hour).  Figure  33  indicates  that  carbon 
particles  were  clearly  separated  from  copper  in  solidified  naphthalene. 

The  vertical  Bridgman  growth  method  was  also  used  to  separate 
mixtures  of  particles.  Carbon  and  copper  particles  were  ul trasonical ly 
dispersed  into  molten  naphthalene  (BNM-3).  The  growth  tube  (10.5  mm 
I.U.)  was  evacuated  and  sealed  at  0.05  Torr.  The  growth  rates  were 
increased  from  10  to  40  mm/hour.  The  solidified  crystal  was  about 
7  cm  long  at  the  end  of  a  3-hour  run.  Most  of  the  carbon  was  trapped 
in  the  upper  portion  of  the  crystal,  while  copper  was  incorporated  in 
the  middle  portion.  Although  the  difference  in  their  individual  V 

c 

was  slight  (about  21  mm/hour  for  carbon  and  18  mm/hour  for  copper)  in 
naphthalene,  they  were  still  separated  from  each  other  by  programmed 
solidification.  A  mixture  of  red  iron  oxide  and  carbon  was  not 
separated  appreciably  in  naphthalene  (BNM-1 ) .  However,  the  growth 
tube  was  not  evacuated  and  a  growth  of  bubbles  on  the  freezing  inter¬ 
face  adversely  affected  the  separation. 

Both  growth  methods  have  shown  experimentally  that  a  mixture 
of  particles  can  be  separated  from  one  another  by  programmed 
solidification  of  a  suspension  in  an  appropriate  liquid.  We  call 
our  invention  "particle  chromatography".  Th's  new  separation  technique 
may  be  used  to  identify  particles  from  a  mixture  of  various  kinds  of 
particles,  derived  from  drugs  and  particulate  pollutants,  as  examples. 
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B.  Size  Classification 

i 

Screen  analysis  is  a  simple  and1 accurate  method  of  fractionating 
materials  into  various  sizes.  It  is  widely  used  to  determine  the  size 
distributions  and  to  yield  fractionated  products  which  may  be  studied 

i 

separately  from  the  whole  material.  .However,  the  method  is  applicable 
only  to  about  44  ym  or  325  mesh.  Below  this  size  th&re  are  several 
methods  to  analyze  the  distribution  of  sizes,  hut  it  appears  that 

I 

only  the  method  of  elutriation  may  be  employed  to  obtain  fractionated 

!  ' 

products  in  the  size  range  of  1  to  50  ym  [51].  '  The  basis  of  the  , 
elutriation  method  is  that  an  upward  velopity  of  fluid  supports  only 
particles  smaller  than  a  given  size.  The  smaller  particles  are 

■  .  i 

carried  upward,  while ‘large  particles  settle  in  the  stream.  The  col- 

I 

lected  particles  may  be  further  elutriated  at  different  velocities 
to  make  finer  size  fractions  and  thus  it  is  difficult  to  start  with 

'I  i 

small  amounts  of  particles.  Here  we  have'  a  method  to  make  size 

I 

fractions  for  a  small  quantity  of  particles.  ,  1 

Previous  studies  [5, 6, 7, 9]  indicated  the  dependence  of  the 

i  '  .  1 

critical  velocity,  on  particle  size.  Larger  particles  have  lower 

V  ,  smaller  particles  have  higher  V  .  Therefore,  particles  with 

C  1  ■,  \  \  ,  ^ 

various  sizes  could  be  sieved  in  a  matrix  material  by  gradually  in- 

1 

creasing  the  freezing  rates  during  controlled  solidification.  A 

}  ! 

mixture  of  silver  powder  and  naphthalene  (BNSC-3)  was  solidified  by 
the  Bridgman  growth  method.  The  freezing  rates  were  increased  from 

i  i  i 

6  to  26  mm/hour.  The  solidified  crystal  (5.2  cm)  was  transversely 
sliced  into  four  equal  sections.  The  silver  particles  in  each  section 


109 


were  collected  on  a  Mi  1 1 i pore  filter,  and  were  examined  under  a 
microscope.  The  size  distributions  varied  from  the  top  to  bottom 
sections.  The  micrographs  of  the  top  and  bottom  sections  are  shown 
in  Figure  34.  Small  particles  were  dominant  in  the  top  section,  but 
some  large  particles  were  also  present  in  this  section,  which  might 
have  been  pushed  by  grain  boundary  grooves  and  grain  boundary  triple 
points  [9].  The  bottom  section  contained  many  large  particles  as 
well  as  small  particles  These  small  ones  were  probably  incorporated 
in  the  grooves  between  the  crystal  and  glass  tube  where  Vc  was  found 
to  be  substantially  reduced  in  this  study.  The  method  could  probably 
be  improved  by  growing  a  single  crystal  with  a  slightly  concave 
freezing  interface  so  that  the  effects  of  grain  boundary  and  grooves 
next  to  the  wall  can  be  eliminated. 

I 

C.  Economics 

The  critical  trapping  velocity  will  have  a  major  influence  on 
the  cost  of  making  particle-free  organic  compounds  since  particles  can 
only  be  separated  from  the  matrix  material  at  freezing  rates  below  Vc< 
Because  the  value  of  Vc  is  different  for  each  particle-organic  system 
pair,  the  cost  will  be  strongly  dependent  on  the  system  being 
purified.  An  organic  compound  usually  contains  several  different 
kinds  of  insoluble  particles,  the  identities  of  which  may  or  may  not 
be  known.  This  makes  estimating  the  cost  of  producing  particle-free 
ultrapure  chemicals  very  difficult  unless  one  knows  (a)  the  kinds  of 
particles  and  their  values  of  Vc,  and  (b)  the  methods  of  determining 
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Larger  particles  dominated  in  the  bottom 
section,  200X. 


(B)  Smaller  particles  dominated  in  the  top 
section,  200X. 


Size  classification  of  spherical  Ag 
particles  (less  than  $0 /m)  by  vertical 
Bridgman  growth  of  naphthalene. 


Ill 


the  amount  of  particles  in  organic  compounds. 

Although  carbon  is  a  common  insoluble  impurity  in  solid  organic 
chemicals,  some  other  kinds  of  particles  are  also  frequently  present. 
To  supply  the  first  information  one  has  to  conduct  experiments  to 
determine  the  minimum  value  of  Vc  among  the  various  kinds  of 
particles  present  before  a  particle-free  compound  can  be  produced. 
Prediction  of  Vc  might  become  possible  if  the  foreign  particles  were 
identified.  The  second  problem  of  determining  the  purity  (in  terms 
of  insoluble  impurity)  may  be  qualitatively  solved  either  by  an 
optical  microscope  or  by  a  filtration  method  as  described  in  Chapter 
III.  There  is  no  accurate  and  convenient  method  available  at  the 
present  time  for  quantitative  measurements  of  particles  in  organic 
chemicals.  With  these  limitations,  we  attempted  to  make  particle- 
free  naphthalene.  The  "native  particles"  present  in  naphthalene  were 
difficult  to  remove  by  a  Fisher  zone- refiner  as  mentioned  in 
Chapter  III,  but  most  of  these  particles  were  separated  by  means  of 
a  horizontal  zone-refiner  at  28  mm/hour  with  5  RPM  tube  rotation. 

This  Vc  is  close  to  the  one  for  separation  of  carbon  particles  from 
naphthalene  (see  Table  8). 

With  our  present  knowledge  of  particle  pushing  and  trapping, 
a  crude  cost  estimate  for  producing  particle-free  naphthalene  is  made. 
Because  the  demand  is  not  yet  sufficient  to  warrant  industrial -scale 
production,  a  laboratory-scale  production  is  more  realistic.  Further¬ 
more,  it  is  recommended  that  a  small  group  with  one  engineer  and  one 
technician  be  created  in  a  chemical  company  which  has  established  the 
market  for  ultrapure  organic  chemicals.  With  these  assumptions,  the 
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detailed  economic  calculations  are  shown  in  Appendix  B.  The  total 
capital  investment  is  about  $68,000  and  the  total  product  cost  is 
about  $65,000  for  an  annual  production  of  3CU  Kg  of  particle-free 
ultrapure  naphthalene.  Therefore,  the  product  cost  is  22  cents  per 
gram,  if  a  reagent  grade  chemical  is  used  as  the  feed  material.  This 
cost  will  be  reduced  if  larger  scale  production  is  employed  [52]. 

It  is  interesting  to  compare  the  present  cost  estimate  with 
the  price  of  zone-refined  organic  chemicals  (removal  of  soluble 
impurities).  Over  100  zone-refined  organic  products  are  available 
comnercially.  A  few  examples  are  given  in  Table  12  [53].  The  price 
of  zone  refined  naphthalene  is  abou  $1.42  per  gram  including  profit, 
purification  cost  and  the  material  cost  (about  0.5  cent  per  gram  for 
the  reagent  grade).  One  does  not  know  the  amount  of  profit  made  by 
the  producers  and  suppliers,  but  the  purification  cost  appears  to  be 
a  substantial  amount  and  is  much  higher  than  the  present  estimate. 

Zone  melting  of  naphthalene  [53-55]  is  typically  conducted  in  a 
vertical  zone-refiner  which  normally  requires  slower  growth  rates 
and  numerous  zone  passes  in  order  to  reach  ultrapure  products.  The 
reason  is  that  the  boundary  layer  next  to  the  freezing  interface  is 
much  thicker  in  vertical  zone  melting  than  in  horizontal  zone- 
ref  ini  ng  with  rotation.  Therefore,  the  developed  new  separation 
process  is  economical  to  product  particle-free  ultrapure  organic 
chemi cals . 
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TABLE  12 

PRICE  OF  SELECTED  ZONE-REFINED  ORGANIC  COMPOUNDS 
FOR  PRIMARY  STANDARDS  [53] 


Price 


Compound 

(for  2-1/2  g 

Acetanilide 

$  3.55 

Anthracene 

4.70 

Anthraquinone 

3.55 

Benzoic  Acid 

3.55 

p-Bromoacetanil ide 

3.55 

p-Chloroacetanil  ide 

3.55 

p-Fluorobenzoic  Acid 

12.25 

o-Iodobenzoic  Acid 

4.70 

Naphthalene 

3.55 

Salicylic  Acid 

3.55 

Succinic  Acid 

3.55 

Sulfanilamide 

3.55 

Vanillin 

3.66 

CHAPTER  VI 


INTERPRETATION 


A.  Pushing  Mechanisms 
1 .  Mass  Transport 

Theoretically,  the  steady-state  pushing  of  a  particle  by  an 
advancing  solid-liquid  interface  demands  a  force  to  prevent  incor¬ 
poration  of  the  particle  into  the  growing  crystal,  and  a  supply  of 
liquid  to  the  contact  area  between  the  particle  and  the  freezing 
interface.  The  repulsive  force  must  equal  the  resultant  of  viscous 
drag  [6,8],  surface  energies  [6],  electrostatic  interaction  [6], 
thermal  forces  [8,9],  and  external  forces  [6,8].  Crystallization 
pressure  (Eq.  (2-3))  is  the  maximum  force  that  can  be  exerted  by  a 
growing  crystal  on  the  particle.  On  the  other  hand,  Corren's  phase¬ 
boundary  force  (Eq.  (2-4))  is  the  minimum  force  needed  to  prevent  the 
particle  from  adhering  to  the  crystal.  Neither  of  these  two  limits 
considers  the  transport  of  liquid  which  may  limit  the  rate  of  pushing. 
In  fact,  the  supply  of  material  practically  determines  whether  a 
growing  crystal  traps  or  repels  a  foreign  particle. 

The  viscous  drag  is  apparently  a  major  force  when  a  particle 
is  in  contact  with  a  horizontal  crystallization  front.  If  the  freezing 
interface  is  regarded  as  stationary,  a  crystallization  flow  [39] 
toward  the  interface  is  generated  by  the  crystal  growth.  The  drag 
force  then  results  from  the  flow  of  fluid  to  the  contact  area  between 
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particle  and  interface.  This  force  interacts  with  the  advancing 

interface  through  a  thin  film,  and  produces  a  repulsive  force  on  the 

interface.  In  other  words,  the  crystal  growth  rate  determines  the 

transport  process  which,  in  turn,  produces  a  pushing  force  (Eg.  (2-24)). 

Their  relationship  is  derived  in  Chapter  II.  Equation  (2-37)  describes 

Vc,  below  which  the  particles  are  pushed  by  the  freezing  interface. 

The  calculation  of  Vc  for  the  present  work  is  shown  in  the  next 
section. 


2.  Heat  Transfer 

Theoretical  developments  for  Vc  described  in  Chapter  II  assumed 
that  thermal  conductivities  of  particle,  melt  and  crystal  are  the 
same,  and  thus  that  there  is  no  effect  of  heat  transfer  on  the  inter¬ 
face  shape  behind  the  particle.  However,  the  thermal  conductivity  of 
the  particle  is  in  general  different  from  those  of  melt  and  crystal. 

For  instance.  Figure  35  shows  that  a  particle  (air  bubble)  is  present 
in  front  of  an  ice-water  interface  [56],  When  the  freezing  interface 
approached  the  air  bubble  fixed  under  a  small  plate,  the  interface 
just  under  it  rose  up  (Figure  35-1,  2,  3,  4  and  5).  The  thermal 
insulating  effect  of  the  air  bubble  (kp  <  y  caused  the  water  under 
the  bubble  to  become  colder  than  other  portions  of  water  at  the  same 
level.  This  caused  the  interface  just  under  the  bubble  to  freeze 
faster  than  other  parts  of  interface,  so  that  a  bump  was  forned  on 
the  advancing  interface.  The  air  bubble  escaped  as  soon  as  the  plate 
above  it  was  inclined  (Figure  35-6).  This  suggested  that  a  thin  film 
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of  water  existed  between  the  bubble  and  the  ice  surface. 

The  schematic  diagrams  of  isotherms  and  heat  flows  are  illus¬ 
trated  in  Figure  36  for  the  cases  of  k  <  k  ,  k  =  k  ,  and  k  >  k  . 

P  36  p  36  p  36 

The  presence  of  a  particle  on  a  solid-liquid  interface  with  kp 
different  from  k£  changes  the  isotherms  and  thus  distorts  the  uniform 
heat  flow  which  would  exist  for  kp  =  k^  or  without  the  particle. 
Because  the  equilibrium  interface  normally  approximates  an  isotherm, 
the  interface  shape  behind  the  particle  becomes  convex  for  k  <  ke 

and  concave  for  k  >  k..  Among  these  three  cases,  the  case  with 

P  36 

kp  >  k^  gives  an  effect  similar  to  that  due  to  the  presence  of  a 
heavier  particle.  Both  can  cause  an  indentation  on  the  interface 
below  the  particle  and  reduce  the  mass  transport.  Experimentally, 
copper  particles  gave  a  smaller  \IQ  than  carbon  particles  in  the 
Bridgman  growth  of  naphthalene. 

The  effect  of  particle  conductivity  on  the  interfacial  tempera¬ 
ture  gradient  is  not  exactly  known  [8,9],  but  the  effect  of  inter¬ 
facial  temperature  gradient  on  the  particle  pushing  is  discussed  here. 
The  heat  flux  qz  across  a  particle  through  a  vertical  central  line 
can  be  expressed  by 


q  =  |<  — Sk  AL  =  ^  AX 
qz  Kp  dz  Kp  Az  Kp  2T 


(6-1) 


The  temperature  difference  between  the  top  T^  and  bottom  T^  of  the 
particle  is  then 


AT  -  Tt  -  Tb  -  ^  q 
P 


(6-2) 
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The  above  equation  can  also  be  used  to  explain  the  interface  shapes 

behind  the  particles  as  shown  in  Figure  36.  At  a  constant  steady- 

state  heat  flux,  AT  is  decreased  as  k  increases.  For  the  case  k  > 

P  P 

kji’  ^b  1S  hl9^er  ,an  ^e  temperature  in  the  liquid  at  the  same  level 
and  thus  the  interface  melts  back  (Figure  36C).  The  reverse  is  true 
for  the  case  of  kp  <  k^  and  the  interface  rises  up  (Figure  36A). 

When  the  imposed  temperature  gradient  in  the  melt  is  increased,  i.e., 
qz  is  increased,  AT  increases.  This  means  T^  increases  while  T^ 
remains  approximately  constant.  In  order  for  q2  to  increase,  either 
T^  must  increase  slightly,  the  interfacial  temperature  must  decrease, 
or  the  thickness  of  the  film  must  decrease.  As  a  result,  the  well 
behind  the  particle  becomes  deeper  for  the  case  of  k  >  k  ,  and  the 

r  ^ 

bump  becomes  smaller  for  the  case  of  k  <  k  .  The  former  is  not 

P  *■ 

favorable  for  particle  pushing,  while  the  latter  may  not  affect  Vc> 

This  interpretation  is  based  only  on  the  above  equation,  and  real 
situations  may  not  be  so  simple.  However,  published  data  [9]  indicated 
that  Vc  for  pushing  of  copper  particles  by  the  ice-water  interface 
(k  >  k  )  was  decreased  from  1.0  to  0.47  pm/sec.  for  65  pm  particles, 

p  X, 

from  1.2  to  0.6  pn/sec.  for  40  pm  particles,  and  from  2.4  to  1.7  pm/ 

sec.  for  10  pm  particles  when  the  temperature  gradient  in  the  water 

was  increased  from  1  to  10°C/cm.  As  shown  in  Table  6,  Vc  was  not 

dependent  on  the  temperature  gradient  for  carbon  particles  pushed  by 

freezing  naphthalene  (k  <  kj. 

P  * 

Since  the  heat  carried  away  into  the  growing  crystal  equals 
that  conducted  from  the  melt  plus  the  latent  heat  liberated  at  the 
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freezing  interface  in  a  vertical  Bridgman  growth,  the  heat  balance 
at  the  interface  may  be  expressed  by 

_ks  (ht)  s,i  =  'ki  (bi)h, i  +  v°sAHf  (6“3) 

where  the  subscript  s  denotes  solid,  l  melt,  and  i  at  the  solid- 
liquid  interface.  Both  interfacial  temperature  gradients  change  with 
the  rate  of  solidification.  As  the  growth  rate  increases,  heat 
generated  by  VpsAHf  gradually  dominates  the  conduction  into  the  solid. 
Thus  (dT/dz)s j  increases  while  (dT/dz)^  .  decreases.  Our  experimental 
observations  indicated  that  the  position  of  the  freezing  interface 
moved  downward  in  order  to  dissipate  the  increasing  latent  heat  more 
easily  into  the  coolant.  Because  the  interface  moved  away  from  the 
heater,  (dT/dz)^  ^  decreased  in  addition.  We  define  "rapid 
solidification"  as  the  situation  when  (dT/dz).  .  approaches  zero  or 
heat  conduction  from  the  melt  becomes  negligible  in  Eq.  (6-3).  Under 
this  condition  the  interfacial  gradient  in  the  growing  crystal  can  be 
directly  calculated  by 

(i)s,is-W,'  <«-«> 

Tien's  temperature  distributions  measured  during  the  normal 

freezing  of  naphthalene  [57]  were  re-evaluated  in  light  of  the  above. 

As  shown  in  Table  13,  the  measured  (dT/dz).  .  approached  zero  at  the 

fast  growth  rate  (66  mm/hour)  and  the  measured  (dT/dz)  .  agrees  with 

s ,  1 

the  value  calculated  by  Eq.  (6-4)  for  two  sizes  of  growth  tubes. 

Thus,  for  rapid  solidification  (dT/dz).  .  may  be  directly  estimated 

5  j  1 
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TABLE  13 

INTERFACIAL  TEMPERATURE  GRADIENTS  FOR 
RAPID  SOLIDIFICATION^1 ) 


Temperature  Gradient,  °C/cm 


Growth  Tube 

run .  I .  D . 

Growth  Rate 
mm/hour 

(2) 

Measured' 

(dT/dz)£<i  (dT/dz)s>i 

Calculated  by 
Eq.  (6-4) 
<dT/dz)sJ 

37 

66 

0 

107 

106.2 

19 

66 

6.3 

no 

106.2 

(^Bridgman  growth  of  naphthalene,  86°C  hot  bath  and  0UC  cooling 
bath  for  both  cases  (Tien  [57]). 

(^Graphically  measured  from  Tien's  distributions  b>  present 
author. 
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from  the  latent  heat  and  does  not  depend  on  the, crystal  dimensions 
(this  is  not  true  in  slow  solidification;  see  Appendix  A). 

If  the  critical  trapping  velocity  falls  within  the  criterion 
of  rapid  solidification,  V  may  not  depend  upon  (dT/dz)„  .  because 
the  latter  is  actually  always  small.  The  latent, heat  of  fusion  is 
very  high  for  naphthalene  (35  cal/g)  and  water  (80  cal/g)  and  thus 
one  more  readily  reaches  the  condition  of  rapid  solidification  than 
with  other  transparent  matrix  materials.  Our  results  on  the  carbon- 
naphthalene  system  and  some  of  Cisse  and  Bolling's  results  on  copper- 
water  system  indicated  that  Vc  was  not  significantly  affected  by 

<dT/dzV.i-  1  ,  ■  ’  ■  ' 

1  '  i 

i  1 

3.  Horizontal  Rotation 

—  ,  l  : 

The  particle  pushing  or  trapping  process  in  horizontal  zone¬ 
refining  with  rotation  is  different  from  that  in  the  vertical 

•  ,  1 

Bridgman  growth  method.  The  former,  pushes  particles  by  a  vertical 
freezing  interface,  the  latter  pushes  particles  by  an  upwa.’d  horizontal 
interface.  (The  interface  is  actually  stationary.  The  motion  referred 
to  is  that  of  the  interface  relative  to  the  ingot.)  Particle  motion 
in  the  rotating  zone  eliminates  two  important  factors.  It  prevents 
particles  from  continuous  contact  with  the  advancing  interface,  and 
also  eliminates  gravitational  effects.  Both  of  theke  are  encountered 
in  the  vertical  Bridgman  technique  and^educe  VQ.  This  is  particularly 
true  for  copper  particles  which  possess  a  very, high  density  (~9  g/cm3). 

i  i 

Therefore,  Vc  was  increased  -at  least  300%  for  copper  in  naphthalene 


1 
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when  the  horizontal  rotation  method  was  employed. 

Fluid  motion  in  the  rotating  zone  was  not  a  simple  solid-body 
rotation;  it  was  complicated  by  the  presence  of  radial  and  axial 
temperature  gradients,  air  bubbles  accumulating  on  the  top,  and  solid 
particles  moving  in  the  melt  under  the  influence  of  gravitational  and 
centrifugal  forces.  In  addition  to  local  stirring  effects,  micro¬ 
scope  observation  revealed  that  two  large-scale  convective  currents 
circulated  from  the  bottom  to  the  top  in  opposite  directions  and  they 
met  around  the  center  of  the  melting  zone.  The  flow  patterns  are 

sketched  in  Figure  37.  These  convective  currents  were  not  stable  and 
fluctuated  periodically. 

The  strength  and  duration  of  these  currents  seemed  to  depend  on 
the  temperature  gradients  and  the  size  of  the  molten  zone.  In  addition 
to  the  convective  currents,  two  fluid  streams  moved  toward  the  freezing 
interface.  The  first  one  was  a  crystallization  flow  [29],  equal  to 
''(Ps/Pt).  The  second  one  was  an  axial  flow  (or  secondary  flow)  result¬ 
ing  from  the  existence  of  density  gradients  and  centrifugal  forces  [58], 
Owing  to  friction,  the  melt  next  to  the  solid-liquid  interface  was 
carried  by  it  and  then  forced  to  the  tube  wall  by  the  centrifugal 
acceleration.  Thus  the  mass  of  fluid  which  was  driven  outward  by 
centrifugal  forces  was  replaced  by  an  axial  flow,  which  was  proportional 
to  where  a  is  the  angular  velocity.  These  three  fluid  streams 

would  have  interacted  with  particles  in  contact  with  the  interface 

when  the  growth  rate  approached  V 
i  c’ 

As  soon  as  the  particles  were  in  contact  with  the  interface, 


Convective 

current 


Resistance 

heater 


Figure  37  Flow  patterns  in  a  horizontal 
rotating  zone  (side  view). 
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some  of  the  particles  would  remain  there  through  friction  and 
adhesion.  The  centrifugal  force  acting  on  the  particles  was  balanced 
by  the  shearing  stress  and  viscous  drag.  After  the  particles  attached 
to  the  rotating  interface,  the  pushing  and  trapping  mechanisms  would 
then  be  similar  to  those  in  the  vertical  Bridgman  method.  In  other 
words,  transport  is  the  limiting  process.  Without  considering  the 
convective  currents,  the  radial  velocity  u  in  the  boundary  layer 
next  to  the  interface  is  proportional  to  rft  and  the  thickness  of  the 
boundary  layer  6  is  proportional  to  vnTtyy  When  the  rotation  speed 
is  increased,  u  increases  while  6  decreases.  This,  in  turn,  serves 
to  improve  the  mass  transport  in  the  contact  area  between  the  particle 
and  the  interface.  Further,  this  mass  transport  is  expected  to  be 
enhanced  when  the  convective  currents  are  considered  which  become 
stronger  in  a  larger  diameter  tube.  This  explains  why  Vc  was 
increased  at  higher  rotation  speeds  (copper  and  iron  oxide  particles) 
and  in  a  larger  diameter  tube  (23  mm  O.D.). 

B.  Calculations  of  V 
_ c 

In  order  to  compare  with  the  experimental  results,  calculation 
°f  vc  is  made  here  by  employing  the  simple  Eq.  (2-37)  without  con¬ 
sidering  effects  of  gravity  and  particle  roughness.  Some  of  the 
constants  such  as  a,  n>  aQ  and  (f>(a)  needed  for  the  calculation  are 

not  readily  available  from  the  literature.  These  constants  are  thus 
estimated  first. 
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The  direct  measurement  of  the  specific  free  energy  a  of  the 
scl ■> d- 1  iqu f d  interface  was  recently  made  for  the  transparent  materials 
camphene,  succinoni tri le,  water  and  white  phosphorous  by  Jones  and 
Chadwick  [59].  Their  experimental  figures  were  in  good  agreement  with 
the  values  predicted  from  Turnbull's  theoretical  equation,  [60,61] 

's*  *  (l  t0  ?)  AV"V  V/3  (6-5! 

where  AH^  is  the  heat  of  fusion,  N  is  Avagadro's  number,  and  V  is  the 
molecular  volume.  Because  the  calculated  using  the  lower  bound  of 
the  constant  3/8  in  the  equation  gave  better  results  for  the  above 

transparent  materials,  this  constant  is  used  to  estimate  o  of 

s  lb 

naphthalene,  salol,  camphor  and  benzophenone.  The  calculated  results 
are  shown  iri  Table  14,  along  with  Jones  and  Chadwick's  data. 

The  physical  properties  of  the  supercooled  thin  film  between 
the  particle  and  interface  are  not  known.  (Data  for  properties  of 
supercooled  bulk  liquid  are  even  meager.)  The  viscosities  of  the 
supercooled  melts  for  salol  [G2]  and  benzophencne  [63,64]  are  shown  in 
Figure  38.  When  the  melts  are  above  the  melting  temperature,  the 
Plots  of  £n  n  versus  1/T  give  straight  lines  following  the  Arrhenius 
relation.  But  these  lines  become  curved  when  the  melts  are  super¬ 
cooled.  The  viscosity  of  the  supercooled  melt  is  significantly  higher 
than  the  melt  at  the  melting  point.  Since  the  supercoolings  of  freez¬ 
ing  organics  for  the  present  work  were  not  measured,  the  viscosities 
at  their  melting  points  are  used  for  the  calculations  of  V 

c  ’ 

Experimental  observations  indicated  that  the  supercoolings  of 
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naphthalene  and  camphor  were  very  small,  but  the  supercoolings  of 
salol  and  benzophenone  seemed  to  be  quite  large. 

The  molecular  volume  Vm  was  obtained  by  dividing  the  molecular 
weight  by  the  density  of  the  1  icr  i  at  the  melting  point.  The 
molecular  diameter  aQ  was  then  estimated  from  the  molecular  volume. 

A  plot  of  <+> (ct)  against  a  with  6  as  a  parameter  is  illustrated 
in  "igure  39.  For  a  given  value  of  6,  4(a)  reaches  a  maximum  value 
at  a  <  1/3  (see  Eq.  (2-39)).  By  combining  Eq.  (2-26)  with  <j>(a),  the 
separation  distance  d  between  the  particle  and  interface  can  be  cal¬ 
culated  by  the  following  equation. 

d/a  =  JfclsU  (6-6) 

0  a(l-a) 

The  calculated  values  of  d/aQ  are  tabulated  in  Table  15.  The 

variation  of  d/aQ  is  less  than  a  factor  of  two  when  the  a  is  increased 

from  0.05  to  0.3.  Since  a  represents  the  curvature  of  the  freezing 

interface  which  is  generated  due  to  the  interaction  between  the 

particle  and  interface,  it  may  approach  a  value  of  0.2  ~0.3  at  the 

onset  of  particle  trapping.  In  this  practical  range  of  a,  the  values 

of  d/a  are  almost  constant  for  a  thin  film  existing  between  the 
o 

particle  and  advancing  interface.  For  the  present  calculations  of 
V  ,  B  =  1  and  4>(a)  max  =  0.335  will  be  used. 

Thus  far  the  constants  needed  for  prediction  of  \IQ  have  been 
estimated.  Calculated  results  of  Vc  are  also  presented  in  Table  14 
for  particles  of  1  and  10  m  diameter.  Since  the  experimental  results 
of  Vc  (Table  7)  were  determined  for  the  smallest  trapping  velocities, 


TABLE  15 


VARIATION  OF  d/aQ 

WITH  a  AND 

d/aQ 

6 

a 

8  =  1 

8  =  2 

8  =  3 

0.05 

4.0 

5.0 

6.0 

0.10 

3.3 

4.3 

5.3 

0.20 

2.6 

3.6 

4.6 

0.30 

2.2 

3.2 

4.2 
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(or  Vc  for  larger  particles),  they  are,  in  fact,  in  good  agreement 
with  the  calculated  values. 

C.  Formation  of  Air  Bubbles 

The  freezing  interface  can  reject  dissolved  gas,  exactly  in 
the  same  way  as  any  -  er  solute  for  which  the  segregation  coefficient 
ko  is  less  than  unity.  As  the  gas  concentration  increases  at  the 
interface,  the  melt  may  become  sufficiently  supersaturated  to  nucleate 
bubbles.  The  bubbles  formed  at  the  interface  grow  by  diffusion  of  dis¬ 
solved  gas  from  adjacent  supersaturated  areas  of  the  melt.  The  bubble 
nucleation  may  be  either  homogeneous  or  heterogeneous.  Thorough  treat¬ 
ment  of  nucleation  theory  for  vapor  to  liquid,  liquid  to  solid,  and 
solid  to  solid  transformation  have  existed  in  the  literature  for  some 
time.  Only  recently  has  such  a  theory  been  presented  for  the  liquid 
to  vapor  transformation  [65-68].  The  purpose  of  this  section  is  to 
modify  this  theory  and  to  make  it  applicable  to  nucleation  of  bubbles 

from  the  melt  containing  dissolved  gases  during  controlled  solidifica¬ 
tion. 

The  homogeneous  nucleation  rate  of  bubbles  developed  by  Hirth 
and  Pound  [65]  was  estimated  by  multiplying  the  concentration  of 
critical-sized  buobles  n  by  the  frequency  of  adding  molecules  w.  The 
equilibrium  concentration  was  corrected  by  a  non-equilibrium  or 
Zeldovich  factor  Z  because  the  process  actually  occurs  at  steady  state 
rather  than  equilibrium.  The  nucleation  rate  equation  is  of  the  form 

J  =  Zwn  cm  sec 


(6-7) 
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★ 

The  value  of  n  is  estimated  from  the  metastable  equilibrium 
between  two  phases  of  differing  Gibbs  free  energy 

n*  =  nQ  e/p  (-AG*/kT)  (6-8) 

where  nQ  is  the  total  concentration  of  molecules  in  the  liquid  phase 
★ 

and  AG  is  the  free  energy  change  for  formation  of  a  critical  bubble. 

The  Gibbs  free  energy  of  formation  of  a  spherical  bubble  in  a 
supersaturated  melt  under  a  prevailing  hydrostatic  pressure  (the 
schematical  diagram  of  bubble  formation  is  shown  in  Figure  40)  is 

AG°  =  4nr3  AGy/3  +  4iTr2o  +  47rr3Ph/3  (6-9) 

in  which  AGy  is  the  volume  free  energy  change  and  o  is  the  surface 
energy.  Since  AGy  is  always  negative  for  a  supersaturated  melt, 

Eq.  (6-9)  exhibits  a  maximum  for  a  critical-sized  bubble  r*  where 
(3AG°/3r)  =  0  and  AG°  =  AG  .  The  relation  between  AG  and  r  is 
obtained  as 

AG*  =  47io(r*)2/3  (6-10) 

★ 

where  r  =  -2o/(AGy  +  Ph). 

The  condition  for  mechanical  equilibrium  of  a  bubble  applied  to 
a  critical-sized  bubble  may  be  expressed  as 


★ 

in  which  P  is  the  total  pressure  in  the  critical  bubble.  An  alterna- 

★  ★ 
tive  form  of  AG  is  obtained  through  elimination  of  r  ,  by  Eq.  (6-11) 
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Solid-liquid 
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Figure  40  Schematic  diagram  of  bubble  formation 
at  a  freezing  interface. 
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AG  =  16tto3/3  (P*  -  Ph)2  .  (6-12) 

★ 

P  represents  the  total  pressure  in  a  critical  bubble.  It  is 
the  vapor  pressure  if  the  liquid  is  a  single  component.  It  is  the 
pressure  or  dissolved  gas  if  the  vapor  pressure  of  the  melt  is 
negligible.  It  is  the  sum  of  partial  pressure  of  dissolved  gas  p* 
and  partial  pressure  of  the  melt  p*  when  both  the  melt  and  dissolved 
air  have  considerable  partial  pressures,  such  as  in  the  present 
experiments  on  naphthalene.  The  last  situation  is  treated  here. 
Assuming  ideal  gas  behavior, 

P  =  Pi  +  P2  •  (6-13) 

The  usual  Kelvin  equation  [23]  deals  only  with  a  pure  liquid 
and  a  nucleus  of  its  vapor.  The  vapor  pressure  p*  of  the  solvent 
inside  the  critical  bubble  was  related  to  the  properties  of  a  dilute 
solution  surrounding  the  bubble  by  Ward,  et  al.  [67].  The  relation 
is  derived  here  for  dilute  ideal  solution.  The  chemical  potentials 
for  the  solvent  in  the  solution  and  vapor  phases  are  given  as  [69], 


ui  =  "  kTx2 

(6-14) 

★ 

u-|  =  u"(p  »T)  +  kT  in  xj 

(6-15) 

where  the  subscripts  1  and  2  refer  to  the  solvent  and  solute  respec 
tively.  The  single  prime  indicates  the  liquid  phase  and  the  double 
prime  indicates  the  vapor  phase.  yQ  is  the  chemical  potential  of 
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the  pure  component  and  x  is  the  concentration  injnole  fraction. 

The  chemical  potential  of  the  pure  component  is  a  function 
of  pressure  and  temperature.  At  constant  temperature,  the  Gibbs- 

.  •  i 

Duhem  equation  reduces  to  ■  ' 

1 

! 

ndy  =  VdP  .  .  !  ,  (6"16) 

0  ,  ,  ■ 

i  •  i 

If  the  pure  solvent  is  assumed  to  be  incompressible^  the  above  equa¬ 


tion  is  integrated  as 

i 

p' (ph.T)  =  u;(p„,t)  +  vj(ph-p„)  . 

:  ’  i  1  I  = 

Equation  (6-16)  is  also  integrated  in  the  vapor  phase  and  becomes 

*  ' 

,  Po(P*’T)  *  u0(P-fT)  +  kT  ln  , 

where  P  is  the  vapor  pressure  of  the  pure1  solvent  across  a  flat 

00 

surface.  Substituting  Eqs.  (6-17)  and' (6-16)  into  Eq.  (6-14)  and 
(6-15),  respectively,  the  chemical  potentials  are 


(6-17) 


(6-18) 


■'  "  "'.(Pm.T)  +  v£(Ph-P J  kT  xj 


J1  Mov 


and 


ii  =  y"(poo’T) +  kT  ln  r +  kT  £n  xi  ' 

*  k*  no 


(6-19) 

i  i 

(6-20) 


across 
i  s 


Since  Pq(Poo,T)  =  yJJtP^T],  the  condition  for  equilibrium 
the  curved  .interface  (yj  =  y!j),  using  Eqs.  (6-19)  and  (6-20), 


•k 

P  x" 

kT  in  -p-^  =  vz(Ph-Pj  -  kT  x^ 

00 
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or 


h  = 


(6-21) 


The  nucleation  of  bubbles  occurs  at  the  freezing  interface 

where  the  concentration  of  dissolved  gases  has  the  highest  value  Ci 

(or  x'2  in  mole  fraction).  The  partial  pressure  of  gas  p*  in  the 

critical  bubble  can  be  evaluated  from  (or  x£)  by  the  Henry's  law 
★ 

P2  =  CiHi  where  H  is  an  empirical  constant. 

Because  the  segregation  of  a  dissolved  gas  by  the  freezing 
interface  is  analogous  to  the  segregation  of  a  solute,  the  segrega¬ 
tion  coefficient  kQ  is  defined  here  as  the  ratio  of  the  gas  concentra¬ 
tion  in  the  solid  to  that  in  the  liquid  at  equilibrium.  Using  a 
constant  value  of  kQ  in  the  quasi-steady-state  stagnant  film 
approach  [70],  the  gas  concentration  at  the  interface  during  fractional 
solidification  is  found  to  be 


ko  +-(1“k0)  exP(" 


I  6V  Ps\ 

V  D 


(6-22) 


where  Cb  is  the  gas  concentration  in  the  bulk  melt,  and  6  is  the  thick¬ 
ness  of  the  boundary  layer.  The  vapor  pressure  of  dissolved  gas  at 
the, interface  is  then 


P2 


(6-23) 


The  importance  of  the  gas  concentration  in  the  bulk  melt  C,  is 

b 


■»  -I.,. 


I  Ji 

obvious,  since  it  directly  affects  the  gas  concentration  at  the 
freezing  interface.  The  value  of  Cb  depends  on  the  solidification 
conditions.  If  the  growth  tube  containing  a  solidified  organic  is 
evacuated  and  sealed  as  described  in  Chapter  III,  the  maximum  number 
of  air  molecules  remaining  in  the  vacuum  space  is  N  „  =  p  r/kT 
W  ere  Py  is  the  pressure  and  V"  is  the  volume  in  the  vacuum  space. 

As  soon  as  the  organic  is  melted  before  beginning  solidification, 
the  amount  of  air  in  the  vacuum  space  increases  to  =  p  V"/kT 
where  p2  is  the  partial  pressure  of  air.  The  difference  N2  -  Nq2 
results  from  transport  of  dissolved  air  from  the  melt  to  the  gas" 
space.  Assume  CQb  is  the  gas  concentration  in  the  solid  (its  volume 
VQ)  before  melting,  and  Cb  is  the  concentration  of  the  bulk  melt  with 
a  volume  V1,  the  material  balance  on  the  gases  is 


P2V"  p  V" 

h~=  CobVo  -  CbV'  * 


(6-24) 


Substituting  p2  =  HfaCb  into  the  above  equation,  the  gas  concentration 
in  the  bulk  melt  is  expressed  as 


P  V"  +  C  UV_  kT 


V'kT  +  H,  V" 

b 


(6-25) 


using  Eqs.  (6-13),  (6-21)  a„d  (6-23),  the  activation  free  energy 
for  nucleation  of  a  bubble  is 


* 

AG  = 


_ rjb _ 

yo-k0)exP(-  f  2j) 


+  exp(w 


(6-26) 


The  growth  process  of  a  bubble  consists  of  a  molecule  (dis¬ 
solved  gas  or  solvent)  in  the  surface  of  the  bubble  evaporating  into 
the  critical  bubble.  The  frequency  of  occurrence  is  equal  to  the 
frequency  of  condensation  of  a  molecule  from  the  bubble. 


*9 

47Tr  v 

w  =  — - - exp 


otP  (4Trr  ^) 
(2-nmkT) 1/2 


(6-27) 


And  the  Zeldovich  correction  factor  is 


3AG*v 


1 67i3kTr*6 


9  *4 

4tt  r  \T 


(6-28) 


where  s  is  the  surface  area  of  a  molecule  at  the  interface,  v  is 

the  Debye  maximum  frequency  of  vibration  of  a  liquid  molecule,  AG* 

vap 

is  the  activational  free  energy  for  evaporation,  a  is  the  condensa¬ 
tion  coefficient,  and  m  is  the  molecular  mass. 

The  factor  (Zw)  has  been  estimated  for  pure  liquids  [71]. 

Hirth,  et  al.  [66]  suggested  that  growth  step  of  dissolved  gas  was 
the  favored  process  due  to  AG^^  for  desorption  of  gas  being  less  than 
that  for  evaporation  of  the  melt.  We  assume  here  that  the  growth 
frequency  is  the  same  for  the  solute  and  solvent,  and  the  value  of 

m  can  be  estimated  by  the  relation  l/i4n  =  1  / /m '  +  1/vfii  where  m  is 

9  A  g 

the  gas  molecular  mass  and  m^  is  the  solvent  molecular  mass. 

An  explicit  expression  for  the  nucleation  of  gas  bubbles  is 
obtained  by  substituting  the  Eqs.  (6-26),  (6-27)  and  (6-28)  into 
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All  parameters  in  the  above  equation  are,  in  principle,  known  or 
experimentally  measurable  for  a  given  system.  The  surface  tension 
a  used  in  the  above  derivation  may  not  be  a  constant  but  will  be 
affected  by  solute  adsorption  at  the  interface.  Mirth  and  Pound  [65] 
pointed  out  in  their  ebullition  (or  boiling)  studies  that  the  size 
of  critical  bubble  (of  the  order  of  a  micron  in  diameter  and  contain¬ 
ing  millions  of  molecules)  was  sufficiently  large  that  macroscopic 
thermodynamic  quantities  might  be  employed  with  some  degree  of 
confidence. 

The  above  nucleation  equation  combining  with  Eq.  (6-25) 
indicates  that  increase  of  Cb  (decreasing  V'  or  increasing  Pv),  or 
6V/D  is  favorable  for  the  formation  of  bubbles.  Our  experimental 
results  qualitatively  agree  with  the  above  theoretical  development. 

In  other  words,  evacuation  of  the  growth  tube  (decreasing  Py), 
stirring  or  rotating  the  melt  (reducing  6)  and  slow  growth  rate  V 
hindered  bubble  occurrence,  while  bubbles  formed  at  the  end  of  the 


Bridgman  runs  due  to  increase  of  by  decreasing  V'. 

Gas  bubbles  were  also  found  in  single  crystal  growth  of 
paratel luri te  TeC^  [72]  by  the  pulling  method.  Formation  of  bubbles 
was  prevented  by  the  slow  pulling  rate  (about  1.5  mm/hour)  combined 
with  high  rate  of  the  crystal  rotation  (40-50  RPM).  This  is  in 
agreement  with  our  theoretical  predictions.  The  interpretations  made 
in  the  referenced  paper  are  unsound. 

Recently  Swanger  and  Rhines  [73]  made  theoretical  calculations 
for  homogeneous  nucleation  of  gas  bubbles  in  the  blood  and  found  that 
homogeneous  nucleation  is  impossible  because  the  initial  pressure  of 
the  gas  must  exceed  1900  ATM.  Maeno's  experimental  observations  [56] 
indicated  that  air  bubbles  did  not  form  on  a  smooth  ice-water  inter¬ 
face  unless  the  interface  was  scratched  with  a  glass  rod.  Further,  air 
bubbles  could  develop  between  the  particles  and  the  ice-water  inter¬ 
face,  depending  on  the  types  and  surface  conditions  of  particles. 

Thus,  formation  of  air  bubbles  was  practically  the  result  of  hetero¬ 
geneous  nucleation.  Our  results  revealed  that  carbon  particles  could 
act  as  nucleation  sites  for  air  bubble  formation,  but  copper  particles 
had  no  nucleating  ability  for  bubbles. 


CHAPTER  VII 
CONCLUSIONS 


At  least  three  objectives  have  been  accomplished  from  the 
research  described  in  this  dissertation.  (1)  Discovery  oc  some 
important  and  interesting  phenomena  leading  to  improved  understanding 
of  the  fundamental  nature  of  particle  pushing  and  trapping  at  a 
freezing  interface.  (2)  The  development  of  a  new  separation  process 
enabling  economical  removal  of  foreign  particles  from  organic 
chemicals.  (3)  The  invention  of  a  new  method  to  separate  mixtures  of 
different  particles  by  solidification,  opening  up  a  potentially  new 

field  --  particle  chromatography. 

The  experimental  results  of  the  vertical  Bridgman  crystal 
growth  indicate  that  interface  shape,  growth  container,  growth 
orientation,  and  the  presence  of  bubbles  all  significantly  influence 
the  critical  velocity  for  trapping.  These  variables,  plus  previously 
established  parameters  such  as  particle  size  and  shape,  particle 
roughness,  grain  boundaries,  density  and  thermal  conductivities,  make 
the  process  of  particles  pushing  and  trapping  extremely  complicated. 
For  example,  particles  were  always  trapped  over  a  finite  length  of 
freezing  crystal,  i.e.,  over  a  range  of  velocities  (21  to  35  mm/hour). 
Among  the  four  organic  compounds  used,  only  naphthalene  can  push 
particles  readily  due  to  its  low  viscosity.  The  effects  of  particle 
properties  on  Vc  appeared  to  be  insignif icent,  although  thermal 
conductivity,  density  and  roughness  are  possible  factors  for  the  same 
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size  particles.  The  measured  values  of  Vc  agreed  with  those  calcu¬ 
lated  by  Eq.  (2-37).  The  Vc  for  the  carbon-naphthalene  system  was 
reduced  from  21  to  18  mm/hour  by  trapping  between  the  glass  wall 
and  a  convex  interface.  The  formation  of  gas  bubbles  at  the  freezing 
interface  was  hindered  by  evacuating  the  growth  container,  stirring 
of  the  melt,  and  use  of  a  slow  growth  rate.  These  observations  are 
in  good  agreement  with  the  theoretical  treatment  of  bubble  nuclea- 
tion  (Eq.  (6-29))  developed  here.  Carbon  particles  provided 
nucleation  sites  for  bubbles,  but  copper  particles  did  not.  The 
interfacial  temperature  gradient  in  the  melt  had  no  observable  effect 
on  V  .  A  simple  theory  based  on  Eq.  (6-2)  predicts  that  increase  of 
temperature  gradient  might  reduce  V  if  k  >  k0,  and  would  have  no 
effect  on  particle  pushing  if  kp  <  k^.  For  rapid  solidification  the 
gradient  in  the  solid  can  be  simply  estimated  from  Eq.  (6-4)  with 
small  and  negligible  heat  conduction  in  the  melt  next  to  the  inter¬ 
face,  regardless  of  the  crystal  dimension. 

Horizontal  zone-refining  with  rotation  enhanced  mass  transfer 
and  eliminated  the  detrimental  effects  resulting  from  bubbles  and 
high  particle  density,  and  thus  significantly  increased  Vc  over  that 
obtained  with  the  vertical  Bridgman  technique.  The  measured  Vc  was 
increased  300%  for  copper  (pg  =  9  g/cc)  and  50%  for  carbon  (pg  = 

2  g/cc)  particles  in  naphthalene.  With  this  new  process  it  is  possible 
to  make  clear-cut  separation  of  particles.  The  cost  estimate  showed 
that  product  cost  is  only  22  cents  per  gram  for  particle-free 
naphthalene.  Particle  and  fluid  motion  in  the  rotating  zone  was 


found  to  be  complex.  Three  fluid  streams  flow  toward  the  freezing 
interface.  They  are  crystallization  flow,  secondary  flow,  and  con¬ 
vective  current.  The  radial  velocity  in  the  boundary  layer  should 
have  a  direct  effect  on  transport  of  liquid  in  the  contact  region 
between  the  particle  and  interface.  This  velocity  is  theoretically 
increased  by  increasing  rotation  speed  and  convective  current.  The 
present  results  indicated  that  the  Vc  was  increased  at  higher 
rotation  rates  for  copper  and  iron  oxide  particles  and  in  a  larger 
diameter  tube. 

Potential  applications  suggested  by  this  work  were  also 
explored.  Two  applications  were  found  to  be  feasible.  Separation  of  a 
mixture  of  different  particles  (particle  chromatography)  and  classifica¬ 
tion  of  particles  were  experimentally  tested. 

It  is  hoped  that  the  results  of  these  studies  will  contribute 
to  better  understanding  of  particle  pushing  and  trapping  by  growing 
crystals,  that  the  new  separation  process  will  be  realized  commer¬ 
cially  and  that  the  potential  applications  will  be  widely  accepted. 

The  suggested  future  work  is  briefly  outlined  in  the  next  chapter. 


CHAPTER  VIII 
SUGGESTED  FUTURE  WORK 

Based  on  the  results  O',  this  dissertation,  further  work  on 

particle  pushing  and  trapping  might  be  continued  along  the  following 
1 i nes . 


In  addition  to  the  Bridgman  crystal  growth  and  zone  milting, 
the  Czochralski  technique  is  widely  used  for  the  growth  of  single 
crystals.  A  seed  is  attached  at  the  end  of  a  rotating  rod.  The 
crystal  grown  on  the  seed  is  pulled  vertically  from  the  melt.  Since 
the  seed  controls  the  orientation  of  growing  crystal,  the  study  of 
the  effect  of  growth  orientation  on  V£  is  possible.  This  technique 
does  not  need  a  growth  container  which  limits  the  rotation  speed  in 
horizontal  zone-refining  due  to  particles  being  trapped  next  to  the 
tube  wall  at  high  rotation  rates.  Further,  sectioning  of  solidified 

products  without  containers  is  easier  for  the  purpose  of  microscopic 
investigation. 


Like  horizontal 
method  also  eliminated 


zone-refining  with  rotation,  the  Czochralski 
settling  effects  which  produce  continuous 


contact  of  particles  with  the  interface  during  the  process  of  particle 
pushing.  However,  the  former  pushes  particles  by  a  vertical  interface; 
the  latter  will  use  a  horizontal  interface  to  reject  particles.  The 


rotation  speed  is  an 


important  factor  in  Czochralski  crystal  growth. 
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and  can  significantly  affect  the  shape  of  growing  interface  and  flow 
patterns  in  the  melt  [3,72,74,75].  At  slow  rates  of  rotation, 
natural  convection  currents  dominate  in  the  melt  and  the  interface  is 
normally  convex.  The  interface  becomes  more  concave  at  higher 
rotation  rates  owing  to  an  axial  flow  (secondary  flow)  resulting  from 
centrifugal  acceleration.  It  would  be  interesting  to  study  how 
rotation  influences  Vc  in  this  geometry.  A  mixture  of  carbon 
particles  and  salol  appears  good  materials  to  start  investigating 
particle  separation  by  this  potential  process. 

B.  Particle  Chromatography 

As  described  in  Chapter  V,  horizontal  zone-refining  with 
rotation  has  proved  to  be  a  promising  technique  for  separation  of 
mixed  particles.  For  instance,  carbon  particles  were  readily 
separated  from  a  mixture  of  carbon  and  copper  particles.  This  new 
technique  should  be  explored  further  for  a  mixture  containing  several 
different  kinds  of  particles,  especially  for  metal  or  semiconductor 
materials.  Since  quantitative  analyses  of  these  materials  have  been 
well  established,  determination  of  separation  efficiency  is  possible 
by  analyzing  each  trapped  band.  Th°  separation  efficiency  is  expected 
to  depend  on  the  freezing  rate,  the  speed  of  rotation  and  the  matrix 
material  selected. 

C.  Crystal  Growth  from  Solution 

This  dissertation  so  far  has  been  concentrated  on  particle 
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separation  for  crystal  growth  from  the  melt.  However,  crystallization 
from  solution  is  a  very  common  and  Important  separation  technique  in 
the  inorganic  chemical  process  industries.  As  mentioned  in  the 
chapter  on  Literature  Studies",  the  evidence  has  shown  that  certain 
crystals  growing  from  solution  have  a  tendency  to  reject  foreign 
particles.  No  information  is  available  at  the  present  time  or.  what 
kinds  of  crystals  tend  to  reject  foreign  particles  more  than  the  other 
kinds  of  crystals.  It  has  been  known  that  some  inorganic  salts  can 
crystallize  in  different  forms  and  habits  under  varied  crystalliza¬ 
tion  conditions.  Sodium  sulfate  crystallizes  Na2S04*10H20  from  the 
aqueous  solution  by  cooling  it  below  32°C,  while  anhydrous  Na2S04  is 
formed  above  32°C  by  evaporation  of  water.  If  the  solution  is  mixed 
with  carbon  particles,  the  results  will  indicate  whether  Na2S04-10H20 
or  anhydrous  Na.,S04  rejects  more  particles.  A  similar  experiment  can 
be  made  for  Na2C03  solution.  Sodium  carbonate  can  crystallize  in 
several  different  forms  as  Na^-iwy).  Na^^O,  Na^^O,  and 
anhydrous  Na^.  The  form  incorporating  the  least  amount  of  particles 
should  be  crystallized  during  purification. 

urther  study  may  be  carried  out  with  a  complex  solution,  such 
as  the  Great  Salt  Lake  brine.  Several  companies  are  interested  in 
recovery  of  saline  minerals  as  well  as  magnesium  from  the  brine  of  the 
Great  Salt  Lake  by  the  solar  evaporation  processes  [76].  Fractional 
crystallization  by  solar  evaporation  is  under  dynamic  conditions. 

Certain  salts  are  very  difficult  to  crystallize,  although  they  should 
be  formed  according  to  the  phase  equilibrium.  Minerals  such  as  NaCfc, 


epsomite  (MgS04* 7H20 ) ,  schoenite  (K2S04*MgS04-6H20) ,  kainite  (KC£- 
MgS04 * 3H20 )  and  carnal lite  (KC£*MgC£2*6H20)  may  be  consecutively 
crystallized  in  the  individual  solar  ponds.  However,  if  a  con¬ 
siderable  amount  of  S04  is  removed  from  brine  at  the  beginning  of 
solar  evaporation,  a  different  crystallization  path  is  followed. 

Salts  of  MaC£,  MgS04-7ll20,  KC£  and  carnallite  are  dropped  out  in  the 
solar  ponds.  These  solar  salts  are  harvested  and  are  further  proces¬ 
sed  in  the  manufacturing  plants  to  produce  KC£  (potash),  K2S04> 

Na2S04,  Mg,  etc.  However,  a  considerable  amount  of  the  insoluble 
materials  (such  as  sands,  clays,  muds,  etc.)  trapped  in  the  salts  have 
to  be  removed  by  several  processing  steps  including  thickening, 
flotation,  centrifugation,  and  filtration.  These  separation  steps 
increase  the  capital  and  operating  costs,  and  thus  influence  the 
overall  process  economics  and  profitability.  If  one  knows  some 
crystals  (likely  kainite)  trap  more  foreign  particles  than  others, 
these  crystals  should  be  avoided  in  process  development.  An  experimen 
tal  study  is  needed  to  identify  these  crystals. 


NOMENCLATURE 


2 

A  Area ,  cm 

a  Molecular  diameter,  cm 

a  Interatomic  distance,  cm 

o 

B  Weight  of  load,  dyne 

b  Height  of  a  monoatomic  layer,  cm 

3  -3 

C  Concentration,  g/cm  or  cm 

3  -3 

C^,  Saturated  concentration,  g/cm  or  cm 

c  A  constant  in  Eq.  (2-6) 

? 

D  Diffusion  coefficient,  cm  /sec. 

d  Separation  distance  or  thickness  of  a  thin  film,  cm 

E(r)  Distance  between  a  concave  interface  and  its  lowest  point 
as  a  function  of  r,  cm 

F  Force,  dyne 

Fc  Viscous  drag  at  a  plane  sink,  dyne 

F(a)  Viscous  drag  at  a  curved  sink,  dyne 

Fp  Force  needed  for  diffusion,  dyne 

G  Gibbs  free  energy,  erg 

H  Biot  number  (=  hR/k ) 

H  Henry's  constant,  dyne  cm/g 

AHf  Latent  heat  of  fusion,  cal/g  or  erg/mole 

h  Separation  distance  between  the  particle  and  interface  at 

any  value  of  r,  cm 

2 

h  Heat  transfer  coefficient,  cal/sec.  cm  °K 

-3  -1 

J  Nucleation  rate,  cm  sec 
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t 


K 

k 

k 

ko 

L 

m 

m 

N 

N 

n 

* 

n 

no 

P 

ph 

* 

P 

P 

00 

p 

V 

P 

R 

R 

r 

ro 

Ar 


Total  curvature, icm 
Botlzmann  constant,  erg/molecule  °K 
Thermal  conductivity,  cal/sec.  cm  °K 
Equilibrium  segregation  coefficient 
Latent  heat  of  fusion,  erg/atom 

I  I 

i  1 

Mass  of  a  molecule,  g/molecule 

Slope  of  liquidus  line  in  a  phase  diagram, °C  cm3/g 

Number  of  molecules  (or  a  number) 

Avagadro's  number,  molecules/mole 

1  • 1  i  i 

A  constant  in  Eq.  (2-5) 

1 

Concentration  of  critical  bubbles,  cm"3 

I 

*  !  O 

Concentration  of  molecules  in  the  liquid  phase,  cm 

O 

Pressure,  dyne/cm 

Hydrostatic  pressure  at  the  freezing  interface,  dyne/cm2 

Total  pressure  in  a  critical  bubble,  dyne/ cm2 

'  2  '  '  ' 

Saturated  pressure,  dyne/ cm 

Pressure  in  the  evacuated  space,  dyne/ cm  . 

?  1  . 

Partial  pressure,  dyne/cm 

Heat  flux  in  the  z  direction,  cal/sec.  cm2 

1 

Radius  of  a  particle  or  crystal,  cm 
Gas  constant,  erg/°K  mole 

Distance  in  a  radial  direction,  cm 

'  1 

Effective  contact  radius  between  a  particle  and  an  inter¬ 
face,  cm 

i  1 

Mean  displacement,  cm  , 


15G 

I  : 

!  I 


i 


! 
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Entropy,  erg/°K 

3 

Volume  entropy  of  fusion,  erg/°C  cm 
Molar  entropy  of  fusion,  erg/°K  mole 

2 

Surface  area  of  a  molecule  at  the  interface,  cm 
Temperature,  °C  or  °K 

Temperature  at  the  top  of  the  particle,  °C 

Temperature  at  the  bottom  of  the  particle,  °C 

Temperature  depression  due  to  curvature  effect,  °C 

Temperature  change  due  to  interaction  of  a  particle 
and  interface,  °C 

Supercooling  for  growth  of  crystal,  °C 
Flow  velocity  in  the  contact  area,  cm/sec. 

Radial  velocity  in  a  boundary  layer,  cm/sec. 

3 

Volume,  cm 

Crystal  growth  velocity,  cm/sec. 

Critical  trapping  velocity,  cm/sec. 

3 

Molar  volume,  cm  /mole 

3 

volume,  cm  /molecule 
Work,  erg 

A  parameter  (=  r2(l-a)/2R) 

Concentration  in  mole  fraction 
Non-equilibrium  or  Zeldovich  factor 
Distance  in  a  z-axis  direction,  cm 


A  parameter  for  the  shape  of  interface. 

Condensation  coefficient 
A  constant  in  Eq.  (2-25) 

Thickness  of  a  boundary  layer,  cm 
Viscosity,  g/cm  sec. 

Angle 

Mobility,  cm/ sec. 

Chemical  potential,  erg/molecule 

Chemical  potential  for  a  pure  substance,  erg/molecule 

Debye  maximum  frequency,  sec-^ 

3 

Density,  g/cm 

O 

Interfacial  free  energy,  erg/cm 

p 

A  parameter  (=  a(l-a)  (3-£n  a)) 

Growth  frequency,  sec"^ 

Angular  velocity,  radian/sec. 


Subscripts 


Sol  vent 

Solute 

Bulk  liquid 

Interface 

Gas  phase 

Liquid  phase 

Original  condition 


p 


External  body  or  particle 
Solid  phase 
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s 


Superscripts 

Critical  nucleus  conditions 
Liquid  phase 
Vapor  or  gas  phase 
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The  influence  of  live  length  and  radius  of  a  crystal  on  the  interracial  temperature  gradient  was  eslimalcd 
by  use  of  a  one-dimensional  heat  transfer  calculation  It  was  found  that  tnc  temperature  gradient  is  sub- 
llanlially  independent  of  crystal  length  only  when  the  Biol  number  hR  k  is  large  twhete  h  is  Ihe  heal  transfer 
coefficient  front  the  cry stal  suiface,  R  is  crystal  radius,  and  k  is  thermal  conductivity). 


1.  Introduction 


2.  Theory 


The  temperature  gradient  ai  the  solid-liquid  inter- 
face  of  a  solidifying  crysial  is  an  important  parameter. 
It  influences  the  microscopic  shape  siability  of  the 
interface,  the  generation  of  dislocations,  the  macro¬ 
scopic  shape  of  the  interface,  eutectic  structures,  and 
free  convection  in  the  adjacent  melt.  In  directional 
solidification  techniques  such  as  Czochralski  and 
Rridgman  growth  the  geometry  of  the  crystal  changes 
during  growth,  i.e.,  the  crystal  grows  longer.  As  the 
crystal  grows  the  inierfacial  tcinperalure  gradient 
changes  as  well. 

In  this  paper  we  estimate  the  interfacial  temperature 
gradient  and  in  particular  examine  its  dependence  on 
crystal  geometry.  A  cylindrical  cry  slal  is  assumed  which 
loses  heat  to  its  surroundings  which  are  al  constant 
temperature  7'..  The  heat  transfer  coefficient  h  for  the 
cylindrical  surface  is  assumed  constant,  bul  different 
from  that  at  the  cold  end  of  the  crystal  /iL*.  The  prob¬ 
lem  is  solved  for  a  non-growing  crystal,  but  the  error 
caused  by  neglecting  the  heal  cairicd  by  the  moving 
crystal  is  estimated  A  onc-duncnsion.il  approximation 
is  employed,  with  Ihe  gcometiy  and  heat  transfer 
parameters  as  shown  tn  fig  I. 


'The  hem  transfer  coefficient  is  defined  as  Ihe  rale  of  heal  lots 
per  unit  area  from  the  surface  divn.ed  by  the  ten  pcrjturc  differ¬ 
ence  between  the  surface  and  the  surroundings  lambient  pav  for 
convective  heat  transport,  heat  shields  or  cooling  surfaces  for 
heat  transport  by  radiation). 


The  r,on-dimcnsionalizcd  differential  equation  go¬ 
verning  heal  transfer  in  a  moving  thin  rod  is  '•') 

cO  ....  vO 
<'n' 


c’o 

_  -2//DWI-, 


(I) 


where  0  =  (T-T,)  (TrT,),  > )  ~  :ll(,  a  -  VpC,R  k 
(Pcdet  number).  //  =  UR  k  (lliot  number  for  cylindri¬ 
cal  surface),  and  r  «  it'/R  *■  l.!R  Here  7"  is  the  tem¬ 
perature  at  distance  ;  from  the  interface.  T,  is  the 
interface  temperature,  R  is  the  radius  of  the  crystal.  I 
is  the  growth  rate,  p  is  the  crystal  density,  Cp  is  the 
heat  capacity,  k  is  the  thermal  conductivity,  and  I.  is 


Fig.  I.  Geometry  and  parameters  of  problem. 
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the  length  of  the  crystal  The  boundary  conditions  are 
•*  1  •*  0  (the  solid-liquid  inlerface),  0  •  I ;  (2) 

••l"  MR  m  t  (the  cuiu  end  of  the  crystal), 

Win- -h-£o  - -iiLo:  (3) 

»t  t  «  0  (infinitesmally  short  crystal).  0  =  I.  (4) 

When  llL  ■  II  the  thermal  conditions  at  the  cold  end 
of  the  ingot  are  the  same  as  on  the  cylindrical  surface, 
while  llL  «»  0  corresponds  to  an  insulated  crystal  end 
•nd  II L  —  oo  to  a  perfectly  cooled  end.  T(L)  =  J\. 

We  have  been  unable  to  solve  these  equations  ana¬ 
lytically  because  of  the  moving  boundary  condition, 
eq.  (3).  Although  a  numerical  solution  could  be  ob¬ 
tained  this  is  thought  to  be  unprolitable  since  the 
equations  only  approximate  reality  Therefore,  we 
first  examine  the  solution  for  a  =  0.  which  is*. 

0  -  sinh  (2//)*  ( r  -  r,)  -t-  ^  j *  cosh  (2//)'  (r  - »,)J 

(Hl  /  2  \*  -,-t 

-  sinh  (2//)*r  4  f-J  cosh(2//)’r  .  (5) 

l-rom  this  the  dimensionless  intcrfacial  temperature 
gradient  is  found  to  be 


+  sinh  (2//)*rJ 


x  sinh  (2//)*r+ cosh  (2//)‘rj  .  (6) 

Curves  calculated  from  this  equation  are  shown  in 
fig.  2.  It  is  seen  that  the  predicted  temperature  gradient 
is  independent  of  length  only  for  II  «  II L  =2.  As  the 
crystal  length  increases,  the  dimensionless  tempera¬ 
ture  gradient  approaches  -(2//)'.  The  fractional 
deviation  from  this  asymptotic  value  ts 

ft  -  f  1(2//)  *  -  (cO/t»0|](2//)-  * 

-[|«  -(<;)' 

x  [7^ sinh  (2/,)‘i+ (1)  cosh (2//)‘tj  .  (7) 

•  Reed  considered  this  problem  Lor  radiation  imo  zero  tempera¬ 
ture  with  a  pcrfc.il>  cooled  end  I /t,  -  r.  r,  0  and  A  -  acT‘, 
where  «  is  the  emisvisiiy  and  a  is  the  Siefan  Ooluminn  eon' 

turn)'). 


’  —  —  d  *0i  _i 

- : 

L  *01-  X  /  - - 

•  I  it 

oi  £  /  _ -r . 

-  * 
005-  /x  ^  I 

:  /  '  : 

/  —  Kj/H  •! 

«“r/  - VH-o 

V  - V"'- 

ca-d - 1 - 1  _  J - u  .  »  1  I 

6  t  l  I  L/Ra  s  t  ? 

ri*.  2.  Fred  ivied  dependence  of  ihe  dimensionless  inicrfacial 
lemperaiure  gradicnl  on  inc  length  of  ihe  cr>sial  and  on  the 
Biol  number  H  liRIK  for  Ihe  cylindrical  surface  for  lfL  n  p 
(cold  end  Insulaledi,  //t  //,  and  //L  x  (pcrfcclly  cooled 

end,  TIL)  —  T,).  eq.  16). 


- 

- *-o. 


a - L-l-LLdUJ  __l_Udiiiq — l_LJ.Illiip*Lrf  t*2a=  J 

001  01  ..  10  IQO 

Mt-hpi 

Fit  3.  Length  of  the  crystal  for  which  ihe  iiiierfaual  lempc^a- 
lure  gradient  deviates  by  1%  from  that  for  a  long  crystal,  for 
II Jll  t  and  0  or  x  (same  rcsull  for  0  and  x |;  eq  Id). 
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fig.  4.  Percentage  error  in  the  interracial  temper ature  gradient 
of  a  tong  cr) stul  caused  by  neglecting  the  heat  carried  by  th¬ 
ermal  motion  due  to  growth;  eq  til) 


From  this  wc  may  find  the  length  of  crystal  for  which 
the  interracial  temperature  gradient  differs  from  the 
asymptotic  value  by  any  chosen  amount.  If  /L  <  I, 
we  note  that 

sinh  (2//)'t  s  coslt  (2//)'t  *  }  esp  [C//)'r3 . 

and  eve  find  that  the  length  at  which  the  error  is /L  is 
given  by 

r,  -  (LIU), 


This  is  plotted  in  fig.  3  for  a  I  deviation  from  the 
asymptotic  gradient,  i.e..  /L  =0.01  rf  may  be  re¬ 
garded  as  the  dimensionless  length  below  which  the 
interfacial  temperature  gradient  depends  on  length. 

The  influence  of  crystal  motion  on  the  interfaci.il 
temperature  gradient  can  only  bo  determined  analyti¬ 
cally  for  a  very  long  crystal  (/.  —  cc).  The  steady  state 
solution  (tfl/ct  =  0)  of  eq.  (I)  for  L  -•  oo  isJ) 

0  ■=  e*p{l[*-(*J  +  K»),J'l).  (9) 

which  yields  the  dimensionless  interfaeial  temperature 
gradient 

(d0/dri.~i[x-(*1  +  8ff)<].  (10) 

Thus  the  fractional  error  in  the  interfaeial  temperature 


gradient  of  a  long  crystal  caused  by  neglecting  the 
motion  of  the  crystal  is  found  to  be 

i  -(8//)* 4  r//)*  ; 

^  |  //)'  |* 
which  is  plotted  in  fig  4 
3.  Discussion 

These  results  show  that  when  a  crystal's  length  ex¬ 
ceeds  its  diameter,  its  interfaeial  temperature  gradient 
depends  on  its  length  only  for  small  values  ( £  I)  of 
the  Biot  number  II  This  is  reasonable  since  the  lliot 
number  is  the  ratio  of  the  surface  to  the  longitudinal 
thermal  conductivity,  or  the  ratio  of  the  case  of  losing 
heat  from  the  surface  to  the  ease  of  conducting  heat 
down  the  crystal  Thus,  for  example  if  a  crystal's 
length  is  twice  its  diameter,  then  for  large  II  ( 2:  0  15) 
its  intctf.vdn<)  temperature  gradient  is  influenced  only 
by  the  thermal  conditions  at  the  cylindrical  surface, 
while  for  small  //(<  0  15)  ’.lie  thernal  cor.dittrr*:  at 
the  other  end  of  the  crystal  al>o  influence  the  inter- 
facial  temperature  gradient.  Values  ol  H  can  be  es¬ 
timated  for  typical  crystal  growth  situations  Foi  good 


Tabu  I 


Comparison  of  the  prevent  results  (using  //t///  -  x,  T, 
H  w*  t  cm,  and  h  **  at  T",*)  wnh  ihow  of  ReeJ1 

-OK. 

t 

Material 

r.c 

Cr.O, 

W 

Melting  point,  /*»  l  K) 

1210 

2531 

3640 

Thermal  eor  Juctiviiy,  k 
(W/cm‘K) 

04 

004 

1.4 

Emisiivity,  e 

02 

OS 

0.36 

Heat  transfer  coefficient,  h 

0  002 

0075 

0  099 

lW  cm*  K) 

Biot  number.  II  ■»  hR.'k 

0  005 

1  87 

0  071 

Present  estimate  of  emu  1 

26 

t  4 

7.2 

length  l.t  be>ond  which 
mtcrfjcul  gradient  con¬ 
stant  tfig  3t  tern) 

Reed's  cntie.il  length,  /.,  tern) 

15 

0  8 

39 

Deviation  of  the  in\rfaeij| 

It) 

9  1 

10  8 

gradient  front  value  for 
»cmi*infinitc  costal  using 
Reed  s  L,  in  eq,  tSl.  U  I".) 
Present  estimate  oT  interfaeial 

121 

4880 

1365 

gradient  for  scnn-infmitc 
crystal.  td  T  d;l<> 
-t2>/)*r,/Kt  K.cntl 

Reed's  estimate  uf  td  T! d.-)© 

76 

3000 

865 

(  K/cm) 

Experimental  IdPdrln  from 

1 1 5 

_ 

- 

ref.  6  t  K/cml 
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conductors  wc  generally  find')  H  <  0.1  for  an 
organic  compound  with  h  -  10  1  W  cm1  C  for  free 
convection  to  air’).  k  10  ’  W,cm  'C.  and  K  «  4 
cm.  we  find  II  *»  4 

At  this  point  il  is  worthwhile  lo  examine  the  validity 
of  the  various  assumptions  ih.it  have  been  made.  It  lias 
previously  been  shown’’'4)  that  the  onc  dinieiiMon.il 
approximation  used  here  is  valid  for  II  <  0  2.  IV ■■ 
tunaicly  it  is  for  these  low  values  of  II  that  our  results 
•re  significant  in  that  a  dependence  of  interfacial  tem¬ 
perature  gradient  on  length  is  predicted  then.  I  or 
H  £0.2  we  expect  our  results  io  be  inaccurate  pri¬ 
marily  in  the  sense  that  the  ntterfacial  tentp.-rai.ire 
gradient  will  be  a  function  of  radial  position  Our  con¬ 
clusion  that  the  inlcrfacial  thermal  conditions  are 
substantially  independent  of  crystal  length  is  expected 
to  remain  valid. 

We  have  assumed  h  is  constant,  whereas  in  almost 
all  real  crystal  growth  systems  l<  will  be  a  function  of 
:  and  I,  both  for  convective  heat  transfer  and  for 
radiative  heat  transfer  When  one  estimates  /.-.  n  is 
done  for  conditions  near  the  interface,  which  have 
the  primary  influence  on  interfacial  tcmperaiuie  gra¬ 
dients  Il  is  of  intcresi  io  compare  the  present  results 
with  those  of  Kecd.  as  shown  in  table  I. 


The  error  caused  by  neglecting  the  heat  earned  by 
the  movement  of  a  growing  crystal  is  estimated  in  cq 
(II).  Ihc  worst  possible  case  is  that  of  a  rapidly 
growing,  insulator  Thus  for  an  organic  with  t'  5 
cm  hr.  /?  u  2  cm.  //  «  2  and  j  ■  1.3.  the  error  is 
about  20*  .  ie.  the  interfacial  tcnipcratuic  gradient 
of  a  long  oiganic  crystal  is  20”.,  lower  than  the  quasi- 
steady  state  estimate. 
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APPENDIX  B 


ECONOMIC  CALCULATIONS 

Following  are  the  detailed  calculations  for  the  cost  estimate 
for  producing  particle-free  naphthalene.  The  assumptions  made  for 
the  calculation  are  as  follows: 

(1)  the  production  group  consists  of  one  engineer  and  one 
technician.  This  group  works  Independently  in  a  large  company. 

(2)  The  engineer  is  responsible  for  research,  development, 
and  administration.  The  technician  is  responsible  for  production 
and  maintenance. 

(3)  The  feed  material  is  a  reagent  grade  chemical,  such  as 

I 

Baker's  analyzed  grade  naphthalene. 

A.  Annual  Proauction 

In  the  two-man  production  group,  the  technician  can  handle 
five  of  six  units  of  horizontal  zone-refiners  with  rotation  in  routine 
production.  Each  unit  is  equipped  with  automatic  control,  not 
requiring  much  attention.  The  growth  tube  is  25  mm  O.D.  and  40  cm  long 
(The  purified  ingot  is  22  mm  in  diameter  and  about  30  cm  long.)  The 
freezing  rate  is  assumed  to  be  25  mm/hour  (  <28  mm/hour).  Because 

I 

two  heaters  are  on  each  refiner,  two  zone  passes  are  assumed  to  ensure 
particle-free  naphthalene. 

I 

Trave1  ktae  “  =  16  hours 
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Zone- ref4 n in g  operation  can  be  conducted  overnight  so  that  most  of 
the  technician's  time  is  spent  in  preparations  for  zone-refining  and 
handling  of  zone-refined  products.  Thus,  daily  production  of  this 
group  is  about  five  purified  ingots.  Assuming  250  working  days  per 
year,  the  annual  production  is 

5  x  (ti  x  2.2  x  30  cm3)  x  (1.145  gr/cm3)  x  250 
=  5  x  (207) (1.1 45)  x  250  a  300  Kg/year  . 

B.  Total  Capital  Investment 

(1 )  Fixed-Capital  Investment 

(a)  Six  units  of  horizontal  zone-refiners  with 
rotation.  Design  and  fabrication  by  the 
production  group  or  modification  of  the 


Fisher  zone-refiner  ($800/unit):  $12,000 

(b)  Building  with  auxiliary  facilities:  30,000 

(c)  Laboratory  equipment:  10,000 

(d)  Indirect  costs  (design,  engineering  and 

contingency):  10,000 

Fixed-Capital  Cost:  $62,000 

(2)  Working  Capital,  10%  of  the  Fixed-Capital 

cost:  $  6,000 

Total  Capital  Investment  (1)  and  (2):  $68,000 


C .  Total  Product  Cost 

(1)  Wages  with  100%  overhead  (for  one  engineer 

and  one  technician):  $50,000 

(2)  Depreciation  (10  years)  $62,000  x  0.1:  $  6,000 


(3) 

Feed  materials  (400  Kg): 

$  3,000 

(4) 

Utilities  (electricity,  water  and  fuel): 

$  1,000 

(5) 

Maintenance  and  repairs: 

$  1,000 

(6) 

General  expenses  (administrative  costs 
research  and  development): 

$  2,000 

(7) 

etcC)- 1an60us  (1nterest.  insurance,  taxes, 

$  2,000 

Total 

Product  Cost: 

$65,000 

Cost  of  Production 

$65,000/300  Kg  =  2 2<t/g. 


